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Abstract  of  Dissertation  Presented  to  the  Graduate  School 
of  the  University  of  Florida  in  Partial  Fulfillment  of  the 
Requirements  for  the  Degree  of  Doctor  of  Philosophy 

PARTIAL  NITRIFICATION  OF  HIGH  NITROGEN 
WASTEWATER  IN  THE 
CONTACT  STABILIZATION  PROCESS 

By 

Byonghi  Lee 
August  1992 

Chair:  Ben  Koopman 

Major  Department:  Environmental  Engineering  Sciences 

The  purpose  of  this  research  was  to  evaluate  the  effect 
of  solids  retention  time  (SRT)  and  sludge  recycle  rate  on 
nitrification  of  high-nitrogen  wastewater  in  the  contact 
stabilization  process,  and  to  develop  and  verify  a computer 
model  that  depicts  the  transformation  of  nitrite  as  well  as 
other  nitrogen  forms  in  wastewater  treatment  processes. 

A bench-scale  contact  stabilization  process  was  operated 
at  a temperature  of  25 °C  with  a feed  of  synthetic  wastewater 
containing  230  g/m3  ammonia  nitrogen.  SRT,  over  a range  of 
4.8  to  8 days,  had  little  effect  on  the  proportion  of  effluent 
nitrogen  which  was  present  in  an  oxidized  form,  but  had  a 
significant  effect  on  the  relative  concentrations  of  nitrite 

vii 


and  nitrate  nitrogen.  Higher  nitrite  concentrations  were 
observed  at  low  SRT.  Sludge  recycle  rate  had  a large  effect 
on  the  proportion  of  effluent  nitrogen  in  oxidized  forms  and 
substantially  influenced  the  relative  concentrations  of 
nitrite  and  nitrate  at  low  SRT.  The  proportion  of  oxidized 
nitrogen  in  process  effluent  varied  from  10%  to  65%  of  total 
nitrogen  when  the  sludge  recycle  rate  was  controlled  between 
20%  and  140%  of  influent  flow  rate.  At  an  SRT  of  4.8  days, 
nitrite  varied  from  12%  to  5%  of  the  total  oxidized  nitrogen 
as  the  sludge  recycle  ratio  was  controlled  between  50%  and 
110%.  At  higher  SRTs,  nitrite  was  low  (3%  or  less  of  total 
nitrogen)  and  did  not  vary  with  sludge  recycle  ratio. 

A computer  model  for  the  contact  stabilization  process 
was  developed  using  rate  expressions  from  the  Activated  Sludge 
Model  No.  1 proposed  by  the  International  Association  on  Water 
Pollution  Research  and  Control  (IAWPRC) . New  rate  expressions 
for  reactions  involving  nitrite  were  added  to  the  IAWPRC 
model.  Ammonia  and  nitrite  half-saturation  coefficients  and 
effective  maximum  specific  growth  rates  for  ammonia  oxidizing 
and  nitrite  oxidizing  biomass  in  synthetic  wastewater  were 
determined  experimentally,  while  values  for  the  other 
biokinetic  and  stoichiometric  model  parameters  were  taken  from 
the  literature.  The  computer  model  gave  predictions  for  total 
oxidized  nitrogen  that  were  very  close  (±  6%  on  the  average) 
to  experimental  results,  but  tended  to  underpredict  nitrite 
concentrations . 
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CHAPTER  1 


INTRODUCTION 

Nitrogen  removal  from  municipal  and  industrial 
wastewaters  is  becoming  increasingly  important.  The  most 
cost-effective  means  of  nitrogen  removal  usually  consists  of 
biological  nitrification  (the  oxidation  of  ammonia  to  nitrate) 
followed  by  biological  denitrification  (the  reduction  of 
nitrate  to  gaseous  end  products) . Frequently,  activated 
sludge  processes  are  applied  because  they  can  be  configured  to 
utilize  the  organic  compounds  present  in  influent  wastewater 
as  electron  donors  for  denitrification. 

Invariably,  activated  sludge  systems  utilized  for 
nitrogen  removal  are  designed  to  achieve  "complete 
nitrification"  of  influent  wastewater.  This  is  true  even  when 
only  partial  nitrogen  removal  is  required,  much  of  the  reason 
for  this  practice  lies  in  the  biokinetics  of  nitrifier  growth, 
which  in  most  activated  sludge  process  configurations,  makes 
nitrification  on  "all  or  none"  phenomenon.  Operation  is 
stable  in  a non-nitrifying  mode  or  in  a fully  nitrifying  mode. 
In  between,  problems  are  experienced  due  to  the  accumulation 
of  a nitrification  intermediate — nitrite — which  leads  to 
erratic  chlorine  demand,  rising  sludge,  and  effluent  toxicity. 
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Gujer  and  Jenkins  (1975)  were  the  first  to  suggest  that 
the  contact  stabilization  process  could  be  used  to  achieve 
partial  nitrification.  They  developed  a mathematical  model 
for  the  process  that  showed  the  rate  of  sludge  recycle  to  have 
a large  effect  on  the  extent  of  nitrification.  Alexander  et 
al.  (1980)  confirmed  the  capacity  of  the  contact  stabilization 
process  for  partial  nitrification. 

Since  nitrification  is  an  energy-intensive  process  and  is 
also  a significant  determinant  of  reactor  volume  requirements 
in  activated  sludge  processes,  it  is,  in  fact,  wasteful  to 
design  for  complete  nitrification  when  only  partial  nitrogen 
removal  is  required.  The  contact  stabilization  process  would 
appear  to  be  a promising  alternative  to  conventional  process 
configurations  in  such  circumstances.  A major  question 
remaining  is  the  propensity  for  nitrite  accumulation,  however. 
None  of  the  previous  research  has  addressed  this  question. 

The  purposes  of  the  present  research  were  as  follow: 

1.  To  determine  the  half-saturation  and  effective 
maximum  specific  growth  rate  coefficients  for 
ammonia  oxidizers  and  nitrite  oxidizers  growing  on 
a high-nitrogen,  synthetic  wastewater. 

2.  To  operate  a bench-scale  contact  stabilization 
process  to  determine  the  effect  of  SRT  and  sludge 
recycle  ratio  on  nitrification. 
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3 . To  evaluate  a computer  model  that  includes  nitrite 
as  a wastewater  component  by  comparing  it  with 
experimental  data. 


CHAPTER  2 
LITERATURE  REVIEW 

Microbiology 

There  are  several  species  of  microorganisms  causing  the 
oxidation  of  ammonia  nitrogen  to  nitrite  nitrogen.  Principal 
organisms  are  the  genera  Nitrosomonas  (N.  europaea  and  N. 
mpnocella)  and  Nitrosococcus . Nitrosolobulus  multiformis  and 
Nitrosospira  briensis  isolated  from  soils  are  also  known  to 
carry  out  the  oxidation  of  ammonia  nitrogen  to  nitrite.  Since 
Nitrosomonas . especially  N.  europaea . isolated  from  wastewater 
treatment  environments  has  a major  role,  the  generic  term 
Nitrosomonas  is  used  to  refer  to  all  ammonia-oxidation 
bacteria.  These  bacteria  are  chemoautotrophic  bacteria  which 
obtain  energy  through  a chemical  oxidation  step  and  use 
inorganic  carbon  such  as  carbon  dioxide,  carbonic  acid, 
bicarbonate,  and  carbonate  as  a carbon  source. 

The  genera  Nitrobacter  (N.  agilis  and  N.  winogradskv)  and 
Nitrosocvstis  are  known  as  the  nitrite-oxidizers  in  wastewater 
environments.  Although  Nitrosococcus  mobilis  and  Nitrosospira 
gracilis  have  been  isolated  from  marine  environments,  nitrite 
oxidizers  are  commonly  referred  to  as  Nitrobacter  (Barnes  & 
Bliss,  1983) . 
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Stoichiometry 

Ammonia  oxidation  has  several  steps  making  hydroxylamine 
and  other  unstable  intermediates.  The  first  step  is  the 
formation  of  hydroxylamine  releasing  little  or  no  energy  as  in 
equation  (2-1) . 

NH4+  + 0.502  > NH2OH  + H+  (2-1) 

The  subsequent  step,  oxidation  of  hydroxylamine  to  nitrite, 
is  known  as  a semi-cyclic  process  involving  the  cytochrome 
system.  A reaction  has  been  summarized  as  in  equation  (2- 
2)  : 


NH2OH  + 02  > H+  + N02'  + H20  + Energy  (2-2) 

This  step  produces  the  majority  of  energy  in  the  oxidation 
of  ammonia.  The  complete  equation  for  ammonia  oxidation  is 
equation  (2-3) . 

NH4+  + 1.502  > N02'  + H20  + 2H+  + (240  - 350  KJ)  (2-3) 

There  are  no  intermediate  products  in  nitrite  oxidation 
which  is  expressed  in  equation  (2-4) 


N02'  + 0 . 502 


> N03'  + (65  - 90  KJ) 


(2-4) 
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If  the  components  of  Nitrosomonas  and  Nitrobacter  are  assumed 
to  be  represented  as  C2H7N02,  the  stoichiometric  equation  for 
cell  growth  is: 


15C02  + 13NH4+  > 10NO2'  + 3C5H7N02  + 23H+  + 4H20  (2-5) 

Nitrosomonas 


5C02  + NH4+  + 10N02'  + 2H20  — > 10N03'  + C5H7N02  + H+  (2-6) 

Nitrobacter 

Carbon  dioxide  is  a carbon  source  for  nitrifying 
bacteria  and,  also,  has  the  buffering  capacity  against  the 
change  of  proton  concentration  in  water  through  changing 
chemical  forms,  e.g.,  carbon  dioxide,  bicarbonate  and 
carbonate.  The  hydrogen  ions  produced  in  equations  (2.1), 
(2,5)  and  (2,6)  can  combine  with  these  carbon  materials, 
therefore,  three  comprehensive  equations  can  be  made: 


NH4+  + 1.502  + 2HC03'  > 

N02'  + 2H2C03  + H20  + (240-350  KJ)  (2-7) 


13NH4+  + 2 3HC03'  > 

8H2C03  + 10NO2'  + 3C5H7N02  + 19H20  (2-8) 


NH4+  + 10NO2‘  + 4H2C03  + HC03'  > 

iono3’  + 3h2o  + c5h7no2 


(2-9) 
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Assuming  the  Nitrosomonas  cell  yield  is  0.15  g VSS/(g  NH4+  - 
N) , equations  (2-7)  and  (2-8)  can  be  combined: 

55NH4+  + 7 602  + 109HC03‘  > 

C5H7N02  + 54N02'  + 57H20  + 104H2C03  (2-10) 

Nitrosomonas 

and  assuming  the  Nitrobacter  cell  yield  is  0.02  g VSS/(g  N02 
- N) , equations  (2-2)  and  (2-9)  can  be  combined: 

400N02‘  + NH4+  + 4H2C03  + HC03‘  + 19502  > 

C5H,N02  + 3H20  + 400N03'  (2-11) 

Nitropacter 

and  the  complete  reaction  for  nitrification  can  be  made  by 
combining  equations  (2-10)  and  (2-11) : 

NH4+  + 1.8302  + 1.98HC03'  > 

0.021C5H7N02  + 1.041H20  + 0.98N03'  + 1.88H2C03  (2-12) 

The  Maximum  Specific  Growth  Rates  and  Half-Saturation 
Coefficients  for  Nitrosomonas  and  Nitrobacter 

The  coefficients  such  as  maximum  specific  growth  rates 
and  half-saturation  coefficients  for  Nitrosomonas  and 
Nitrobacter.  respectively,  have  important  roles  in  a model  for 


nitrification. 
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A pure  Nitrobacter  culture  was  studied  to  evaluate 
several  kinetic  parameters  in  nitrification  (Gould  and  Lees, 
1960) . Williamson  and  McCarty  (1975)  proposed  the  infinite- 
dilution  method  to  estimate  half-saturation  coefficients  for 
Nitrobacter  and  Nitrosomonas  using  pure  cultures  of  both 
microorganisms.  The  Lineweaver-Burk  procedure  was  applied  to 
estimate  half-saturation  coefficients  in  their  study.  A model 
for  contact  stabilization  was  adopted  to  estimate  the  maximum 
specific  growth  rate  for  Nitrosomonas  (Gujer  and  Jenkins, 
1975) . The  comparisons  of  experimental  and  simulation  results 
gave  the  maximum  specific  growth  rate  of  Nitrosomonas 
(Alexander  et  al.,  1980).  Eric  and  Murphy  (1980)  showed  the 
different  maximum  specific  growth  rates  for  Nitrobacter  and 
Nitrosomonas  depending  on  different  influent  Total  Kjeldahl 
Nitrogen  (TKN)  concentrations  and  used  a first  order 
differential  equation  to  compare  the  evaluated  values.  2- 
chloro-6-trichloromethyl  pyridine  (N-Serve)  was  used  to 
inhibit  Nitrosomonas  activity  when  the  specific  growth  rate  of 
Nitrobacter  was  evaluated.  Temperatures  were  not  specified  in 
their  study.  The  infinite  dilution  method  (Williamson  and 
McCarty,  1975)  was  tested  and  compared  with  experimental 
results  using  a respirometer  by  Cech  et  al.  (1985)  . They 
concluded  that  the  infinite  dilution  method  had  more  accuracy 
than  their  method  in  the  estimate  of  the  half-saturation 
coefficient  for  Nitrosomonas . Philbrook  and  Grady  (1985) 
proposed  the  modified  infinite-dilution  method  to  estimate  the 
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half-saturation  coefficient  for  the  degradation  of  2- 
chlorophenol . This  method  will  be  presented  in  Chapter  3 in 
this  study. 

Alleman  (1984)  presented  maximum  specific  growth  rates 
and  half-saturation  coefficients  for  Nitrosomonas  and 
Nitrobacter . Terashiman  and  Ishikawa  (1984)  compared 
simulations  and  experimental  results  in  an  oxidation  ditch  and 
presented  the  half  saturation  coefficient  for  nitrification. 
Chudoba  et  al.  (1985)  investigated  the  nitrification  kinetic 
parameters  according  to  reactor  types  and  found  that  the  half- 
saturation coefficient  and  the  maximum  substrate  utilization 
rate  were  higher  in  a plug  flow  (PF)  reactor  than  in  a 
completly  mixed  (CM)  reactor.  The  average  value  for  the  half 
saturation  coefficient  for  nitrification  in  a PF  reactor  was 
0.65  g/m3,  as  compared  to  0.37  g/m3  in  a CM  reactor.  Henze  et 
al.  (1986)  introduced  several  kinetic  parameters  in 
nitrification  along  with  an  activated  sludge  model.  The 
maximum  specific  growth  rate  and  the  half  saturation 
coefficient  for  ammonia  oxidizers  were  0.8  day'1  and  1.0  g NH3- 
N/m3  at  2 0 "C  and  0.3  day'1  and  1.0  g NH3~N/m3  at  10  °C, 
respectively.  Antoniou  (1989)  investigated  the  maximum 
effective  specific  growth  rate  and  the  half  saturation 
coefficient  for  nitrification  in  an  activated  sludge  from 
Kanapaha  water  reclamation  facility,  Gainesville,  FL.  He 
found  the  maximum  effective  specific  growth  rate  had  a maximum 
value  of  0.564  1/d  at  20°C  and  average  pH  of  8.2.  Half 
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saturation  coefficient  was  reported  as  0.2  g NH3-N/m3  at  a 
temperature  of  20 °C  and  a pH  of  7.8.  Gee  (1990b)  found  the 
maximum  specific  growth  rates  and  the  half  saturation 
coefficients  for  Nitrosomonas  and  Nitrobacter . Yoshioka  et 
al.  (1982  a,  b)  evaluated  kinetic  parameters  of  nitrifiers, 
e.g.,  Nitrosomonas  and  Nitrobacter.  isolated  from  mud  in  a 
lake  using  direct  counting  method  in  a pure  culture  at  2 0 "C 
temperature  and  pH  7.711.0.  They  determined  the  maximum 
specific  growth  rates  and  the  half  saturation  coefficients  for 
both  microorganisms  and  found  ammonia  oxidizers  in  a mixed 
culture  had  the  similar  value  of  maximum  specific  growth  rate 
as  in  a pure  culture.  Table  2-1  shows  the  summary  of  both 
kinetic  parameters. 


Kinetics  and  Models  in  Nitrification 


The  nitrification  could  be  divided  into  two  steps: 
ammonia  oxidation  and  nitrite  oxidation.  Equation  (2-13) 
shows  each  reaction  step. 

rNH4  rN02 

NH4+  > N02‘  > N03'  (2-13) 

rNH4  and  rN02  represent  the  oxidation  rate  of  ammonia  and 
nitrite  nitrogen,  respectively.  Though  two  reaction  rates, 
rNH4  and  rN02,  are  considered  as  a zero  order  (Wong-Chong  et 
al.,  1975),  the  Monod  type  model  has  been  widely  accepted  to 
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express  each  reaction  (Sharma  & Ahlert,  1977) . The  expression 
of  a Monod  equation  for  the  growth  of  Nitrosomonas : 


r 


s 


v 

NH4A 


^ANH4  SNH4 


V K 4-  ^ 

XNH4A  nH4A  °NH4 


(2-14) 


where,  rs  = ammonia  oxidation  rate,  g/m3  day 

Am  = specific  growth  rate  of  Nitrosomonas . day'1 

#ANH4  = maximum  specific  growth  rate 
for  Nitrosomonas . day'1 

SNH4  = concentration  of  ammonia  nitrogen,  g/m3 

Ynh4a  = yield  coefficient  for  ammonia 
oxidizing  bacteria,  g/m3 


K„H4A  = half-saturation  coefficient  for  ammonia 
oxidizing  bacteria,  g/m3 


Equation  (2-14)  can  be  used  to  express  the  rate  of  nitrite 
oxidation  in  terms  of  Nitrobacter  with  nitrite  as  a substrate. 

£ANh4'  SNH4'  KNH4A  and  ynh4a  wil1  be  replaced  with  £AN02,  Sno2,  Kno2a 
and  Yno2a,  respectively,  when  nitrite  oxidation  is  considered. 

Since  dissolved  oxygen  is  frequently  present  at  a level 
low  enough  to  limit  ammonia  oxidation,  the  growth  rate 
equation  should  consider  the  effect  of  dissolved  oxygen.  The 
rate  equation  for  ammonia  oxidation  including  the  effects  of 
ammonia  nitrogen  and,  also,  dissolved  oxygen  is: 


rs  = 


-NH4A 


*ANH4 


L NH4A 


nH4 


K, 


NH4A 


NH4 


+ S 


K0NH4  + S0 


(2-15) 
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Table  2-1.  Characteristic  parameters  for  nitrifying  bacteria. 


Process  AanJ 

(day  ')  (mg 

Knh4a1 2  Environment  References 

NH4-N/L) 

Nitrification 

0.5 (20°C) 

P/C3 

(Willamson, 1978) 

(Biomass:  0.65(21°C) 

A/S4 

(Gujer , 1975) 

Nitrosomonas)  0 . 34  ( 12  0 C) 

A/S 

(Gujer,  1975) 

0 . 65 ( 2 0 ° C) 

A/S 

(Dold,  1980) 

0 . 40 (20  ° C) 

2 . 0 (20  0 C) 

A/S 

(Stall,  1974) 

0.33  ( 2 0 0 C) 

1. 0 (20°C) 

A/S 

(Downing,  1964) 

1.08 (20°C) 

0.063  ( 2 0 ° C) 

A/S 

(Poduska,  1975) 

1 . 50 ( 3 0 ° C) 

P/C 

(Engel,  1958) 

0 . 55 ( 20  ° C) 

A/S 

(Alexander,  1980^ 

o 

• 

Ul 

N) 

o 

o 

n 

A/S 

(Beccari,  1979) 

0.96-1.92 

0.6-3 . 6 

A/S 

(Alleman,  1984) 

0. 6 (20°C) 

A/S 

(Warner,  1986) 

0.8 (20”C) 

1 . 0 (20  ° C) 

A/S 

(Henze, 1986) 

0.3 (10°C) 

1. 0 (10°C) 

A/S 

(Henze,  1986) 

0.28 (15°C) 

0.2 (20  °C) 

A/S 

(Antoniou,  1989) 

0 . 56  ( 2 0 0 C) 

A/S 

(Antoniou,  1989) 

0 . 576 ( 23  ° C) 

0.7 (23°C) 

P/C 

(Gee,  1990) 

0.58  ( 2 0 ° C) 

2 . 8 ( 2 0 ° C) 

A/S 

( Yoshioka , 1982 ) 

0.7  (PFS) 

A/S 

(Chudoba,  1985) 

0.37  (CM6) 

A/S 

(Chudoba,  1985) 

2 

.59-4.59 (20°C) 

SRW7 8 

(Stratton,  1967) 

0.65(20°C)  0.6(20°C)  TSWB  (Knoles,  1965) 


1 maximum  specific  growth  rate  for  Nitrosomonas 

2 half-saturation  coefficient  for  Nitrosomonas 

3 pure  culture 

4 activated  sludge 

5 plug  flow  reactor 

6 completed  mixed  reactor 

7 synthetic  river  water 

8 Thames  estuary  water 
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Table  2-1.  continued 


Process 

^ANO?. 

(day1)  (mg 

K 10 

n02A 

N02-N/L)  Environment  References 

Nitratif ication 

0.14 ( 2 0 ° C) 

A/S 

(Downing,  1964) 

(Biomass : 

0. 07 (20°C) 

P/C 

(Willamson,  1978) 

Nitrobacter) 

1.44 (20°C) 

A/S 

(Poduska,  1975) 

1. 39 (30°C) 

22 . 0 (20  °C) 

P/C 

(Boon,  1962) 

(0.34 

-1.77  ( 2 0 0 C) 

SRW 

(Stratton,  1967) 

0.84 (20° ) 

1.9 (20°C) 

TEW 

(Knoles,  1965) 

1.99 (28°C) 

5.01(28°) 

P/C 

(Gould, 1960) 

0.27  ( 2 0 ° C) 

A/S 

(Beccari,  1979) 

0.48-1.44 

0.3-1. 7 

A/S 

(Alleman,  1984) 

0. 14 (20°C) 

A/S 

(Downing,  1964) 

0.07  (20  ° C) 

P/C 

(Willamson, 1978) 

1.44 (20°C) 

A/S 

(Poduska,  1975) 

0.42  ( 2 3 ° C) 

1. 02  (23  °C) 

P/C 

(Gee,  1990  b) 

0 . 8 1 ( 2 0 ° C) 

0. 14 (20°C) 

A/S 

(Yoshioka,  1982) 

1. 39  (30°C) 

22 . 0 (20°C) 

P/C 

(Boon,  1962) 

(0.34 

-1 . 77 , 2 0 ° C) 

SRW 

(Stratton,  1967) 

0.84  (20° ) 

1 . 9 ( 20  ° C) 

TEW 

(Knoles,  1965) 

1.99  (28  ° C) 

5.01(28°) 

P/C 

(Gould, 1960) 

0. 27 (20°C) 

A/S 

(Beccari,  1979) 

0.48-1.44 

0.3-1. 7 

A/S 

(Alleman,  1984) 

0.42  (23  ° C) 

1.02  (23  °C) 

P/C 

(Gee,  1990  b) 

0 . 81  ( 2 0 ° C) 

0. 14 (20°C) 

A/S 

(Yoshioka, 1982) 

9 maximum  specific  growth  rate  for  Nitrobacter 

10  half-saturation  coefficient  for  Nitrobacter 
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Where,  SQ  = concentration  of  dissolved  oxygen,  g/m3 

= oxygen  half-saturation  coefficient  for 
Nitrosomonas . g/m3 


When  nitrite  oxidation  is  considered,  KQN02  will  replace  KqNH4 
and  SN02  will  replace  SNH4.  In  order  to  use  these  kinetic 
equations  to  simulate  both  reactions,  i.e.,  ammonia  and 
nitrite  oxidations,  coefficients  in  equation  (2-15)  should  be 
estimated. 

Based  on  equation  (2-15) , the  model  for  nitrification  was 
developed  (Ekama  and  Marais,  1978;  Haandel  et  al.,  1982; 
Henze,  1986;  Lijklema,  1973;  Water  Research  Commission,  1984; 
Gee  et  al.,  1990).  Lijklema  (1973)  developed  a model 
simulating  the  growth  of  heterotrophic  bacteria,  the  growth  of 
a predator  and  nitrification.  His  work  showed  an  agreement 
between  results  of  experiments  and  simulations.  The  Marais- 
Ekama  activated  sludge  model  (Ekama  and  Marais,  1978)  was 
verified  by  Nichollas  (1982)  and  showed  a good  agreement 
except  the  values  of  TKN  and  MLSS  concentrations  in  an 
aeration  tank.  The  discrepancy  was  explained  as  being  due  to 
the  difficulty  in  the  estimating  nonbiodegradable  particulate 
COD  matters.  Warner  et  al.  (1986)  carried  out  experiments  for 
the  digestion  of  wasting  activated  sludge  using  an  anoxic- 
aerobic  system  and  compared  the  simulations  data  from  the 
model  proposed  by  Alexander  et  al.  (1988)  . Henze  et  al. 
(1986)  developed  a model  combining  carbon  oxidation, 
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nitrification  and  denitrification.  Thirteen  components  in 
wastewater  and  8 reactions  were  incorporated.  The 
characterization  methods  for  the  constituents  in  wastewater 
and  the  procedures  for  estimating  parameters  in  the  model  were 
presented  in  their  work.  Bidstrup  and  Grady  (1988)  wrote  a 
computer  program  for  a model  proposed  by  Henze  et  al.  (1986). 


Nitrite  Build-up  in  Nitrification 
When  the  rate  of  ammonia  oxidation  is  greater  than  that 
of  nitrite  oxidation  (rNH4  > rN02  in  Equation  (2-13)),  nitrite 
starts  accumulating  in  an  oxidation  tank.  Warington,  cited  in 
Anthonisen  et  al.  (1976),  observed  nitrite  in  1876,  before  the 
identification  of  nitrifying  organisms.  When  a soil  column 
was  used  to  oxidize  a higher  concentration  of  ammonia 
nitrogen,  e.g.,  5000  mg  NH4+-N/L,  nitrite  was  accumulated 
(Stojanovic  and  Alexander,  1957) . Free  Ammonia  (FA)  and  Free 
Nitrous  Acids  (FNA)  inhibited  Nitrobacter  activity  and  caused 
the  nitrite  build-up  (Butt  and  Lees,  1960;  Prakasam  and  Loehr, 
1972  ; Wong-Chong  and  Loehr,  1975;  Anthonisen  et  al . , 1976; 

Anthonisen  et  al.,  1977;  Verstrate  et  al.,  1977;  Wong-Chong 
and  Loehr,  1978;  Chudoba  et  al . , 1985).  Anthonisen  et  al. 

(1977)  made  several  diagrams  showing  inhibitory  concentrations 
of  FA  and  FNA  on  ammonia  and  nitrite  oxidations  at  different 
pH.  Verstrate  et  al.  (1977)  found  lower  concentrations  of  FA 
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(0.29  - 0.3  g FA/m3)  caused  the  inhibition  to  the  oxidation  of 
ammonia  than  those  (0.1  - 1 gFA/m3)  reported  by  Anthonisen  et 
al.  (1977).  Fort  et  al.  (1980)  investigated  the  inhibitory 
effects  of  FA  and  FNA  on  nitrification  and  developed  a couple 
of  equations  showing  concentrations  of  FA  and  FNA  as  a 
function  of  temperature  and  pH.  Their  experiment  showed  that 
the  inhibitory  effects  on  nitrification  occurred  at  the  same 
range  of  the  FA  and  FNA  concentrations  reported  by  Anthonisen 
et  al.  (1977).  Alleman  (1984)  investigated  the  nitrite  build- 
up in  nitrification  at  Sequence  Batch  Reactors  (SBR)  and  found 
FA  was  a predominant  inhibitor  under  several  conditions,  e.g., 
a reduced  temperature,  limiting  oxygen  and  carbon  dioxide 
supply,  elevated  pH  and  excess  solids  wasting.  Nitrite 
accumulated  at  a concentration  of  FA  as  low  as  0.05  gFA/m3. 
Randall  et  al.  (1984)  studied  the  effects  of  temperature  on 
nitrite  build-up.  They  found  nitrite  accumulated  when  the 
temperature  was  lower  than  16 "C.  They  concluded  the 
temperature  impact  is  higher  on  Nitrobacter  than  on 
Nitrosomonas ; therefore,  nitrite  can  accumulate  under  low 
temperature  conditions.  Turk  and  Mavinic  (1986)  tried  to 
accumulate  nitrite  in  a series  of  aeration  tanks  with  high 
nitrogen  influent  and  found  initially  accumulated  nitrite  was 
oxidized  after  all  despite  a same  FA  concentration  was 
maintained  in  aeration  tanks.  They  attributed  this  phenomenon 
to  the  acclimation  of  biomass  to  high  concentrations  in 


influent. 
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Nitrification  in  Contact  Stabilization  System 
Contact  stabilization  was  used  to  treat  domestic 
wastewater  based  on  empirical  observations  (Uillrich  and 
Smith,  cited  in  Gujer  and  Jenkins  (1975)).  Gujer  and  Jenkins 
(1975)  adopted  a term  "Nitrification  efficiency"  - the  ratio 
of  nitrate  in  contact  tank  to  that  in  reaeration  tank-  to 
explain  the  partial  nitrification  in  contact  stabilization 
process.  They  assumed  the  nitrate  in  reaeration  tank 
represented  the  total  amount  of  oxidized  nitrogen  by 
autotrophic  biomass.  They  made  a graph  illustrating  the 
nitrification  efficiency  with  respect  to  recycle  ratio. 
According  to  this  graph,  the  nitrification  efficiency 
increases  exponentially  when  sludge  recycle  ratio  is  lower 
than  100%,  whereas  the  efficiency  reaches  the  plateau  when 
sludge  recycle  ratio  is  higher  than  200%.  however,  their 
experiments  did  not  display  the  effect  of  recycle  ratio  on 
nitrification  efficiency  since  the  recycle  ratios  in  process 
operation  were  near  or  above  100%.  And  experimental  data  did 
not  demonstrate  the  effect  of  temperatures  on  nitrification 
efficiency  within  reactor  operation  temperatures,  12  °C  and 
20 "C.  The  change  of  contact  fractional  mass  did  not  show  the 
effects  on  nitrification  efficiency.  Alexander  et  al.  (1980) 
operated  a laboratory  scale  contact  stabilization  process  and 
evaluated  an  activated  sludge  model  proposed  by  Dold  et  al., 
cited  in  Alexander  et  al . (1980).  The  performance  of  a 
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contact  stabilization  process  was  recorded  at  two  different 
SRTs  under  a constant  temperature,  2 0°C,  and  a constant 
recycle  ratio,  200%.  At  6 days  of  SRT,  nitrification 
efficiency  was  75%,  whereas  75.8%  at  10  days  of  SRT. 


CHAPTER  3 


DETERMINATION  OF  AMMONIA  AND  NITRITE 
HALF -SATURATION  COEFFICIENTS  FOR  NITRIFICATION 


Introduction 

An  infinite  dilution  method  was  employed  for  the 
estimation  of  substrate  half-saturation  coefficients  for 
ammonia  (KNH4A)  and  nitrite  oxidizers  (KN02A)  . Williamson  and 
McCarty  (1975)  first  used  this  method  to  measure  the 
coefficient  by  employing  the  Lineweaver-Burk  equation  to 
analyze  their  data.  However,  because  their  technique  of 
analysis  required  three  or  more  trials  to  provide  a reliable 
estimate  of  the  half-saturation  coefficient,  the  modified 
infinite  dilution  procedure  was  used  instead.  Philbrook  and 
Grady  (1985)  used  the  modified  infinite  dilution  procedure  to 
determine  the  half-saturation  coefficient  for  the  degradation 
of  2-chlorophenol  from  the  results  of  two  consecutive 
experiments . 

Development  of  an  Expression  for  Measurement  of 
Kn h4a  and  K|J02a  Using  Fed-Batch  Cultures 

When  a microbial  culture  is  placed  in  a batch  reactor  and 
substrate  is  added  continually  at  a rate  less  than  the  maximum 
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substrate  utilization  rate  at  which  the  culture  can  use  it, 
the  culture  will  regulate  the  substrate  concentration  to  a 
value  consistent  with  the  rate  of  addition.  If  the  volume  of 
feed  added  is  small  with  respect  to  the  reactor  volume  and  if 
the  biomass  formed  over  the  duration  of  the  experiment  is 
small  with  respect  to  the  mass  carrying  out  the  removal,  then 
a pseudo-steady-state  will  be  reached  during  which  the 
substrate  concentration  in  the  reactor  is  constant. 


From  Monod  kinetics, 


U.  X 


B,  A 


= u.  X 


B.A 


where  UA  = specific  substrate  removal  rate 


(3-1) 


= Mamma  / yawha  for  nitritif ication, 

g NH4+1-N/(g  Nitrosomonas  d) 

= ^ 2 / Y , for  nitratification, 

g N02_1-N/(g  Nitrobacter  d) 

Yanh4  = yield  coefficient  for  Nitrosomonas . 

g cells  formed  per  g NH4+1-N  oxidized 

Yan02  = yield  coefficient  for  Nitrobacter. 

g cells  formed  per  g NO^’-N  oxidized 

UA  = maximum  specific  substrate  utilization  rate  of 
either  nitritif ication  or  nitratification 


XB  A = nitrogen  oxidizing  bacteria,  g cells/m3 
= XgANH4  for  Nitrosomonas 
= XBA|g02  for  Nitrobacter 

SN=  concentration  of  nitrogen  components,  g N/m3 
= SNH,  ammonia  nitrogen  for  nitritif  ication, 
g NH4+1-N/m3 

= SH02,  nitrite  nitrogen  for  nitratification, 
g N02'1-N/m3 
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Knh=  half-saturation  coefficients,  g substrate/m3 
= KNH4A  for  Nitrosomonas , g NH,+1-N/m3 
= KN02a  for  Nitrobacter . g N02  -N/m3 


Equation  (3-1)  can  be  rearranged  as  follows: 

XB,A 

Knh  = ( - 1)  SN  (3-2) 

UA  XB,A 

the  substrate  concentration,  SN,  can  be  measured  by  sampling 
from  the  reactor.  Expressions  for  UA  XB  A and  UA  XB  A can  be 
developed  from  a mass  balance  on  a fed  batch  reactor  (see 
Figure  3-1  for  a schematic  of  a fed  batch  reactor) : 


d(V  SN) 

Q SN0  " UAXB,A  V = 0-3) 

dt 

Where  Q = feed  flow  rate,  m3/h 

SN0  = feed  substrate  concentration,  g/m3 
V = bioreactor  volume,  m3 

When  reactor  volume,  V,  is  large  relative  to  the  feed  rate, 
equation  (3-3)  may  be  expressed  as: 


Q S 


NO 


UA  XB.A  V 


dSN 

dt 


(3-4) 


or 
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NOTE:®  = AIR  SCRUBBER 

© = CCfe  GAS  CYLINDER 
® = MIXED  LIQUOR 
© = Cfe  GAS  CYLINDER 


SP)=  SYRINGE  PUMP 
► = GAS  FLOW 


P)  = PRESSURE  REDUCER 
(§)=  ROTOR  METER 
©=  STIRRER 
% = WATER  BATH 

0=  DIFFUSER 

► = LIQUID  FLOW 


Figure  3-1.  A schematic  of  a fed-batch  reactor. 
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Q s 


No 


dS. 


v 


dt 


+ UA  XB,A 


(3-5) 


If  sN  >>  Knh  and  the  substrate  application  rate  (Q1  SN0/V) 
exceeds  the  culture's  utilization  capacity,  then  UA  is  equal 
to  UA.  Rearranging  equation  (3-5) : 


UA  XB,A 


Q,  S 


No 


dS, 


V 


dt 


(3-6) 


If  the  substrate  addition  rate  is  less  than  the  culture's 
utilization  capacity,  a pseudo-steady-state  will  eventually  be 
attained.  In  this  case,  dSN/dt  = 0,  and 


Q,  S 


UA  XB,A 


No 


V 


(3-7) 


On  the  basis  of  these  equations,  it  is  possible  to 
determine  the  for  nitrifying  cultures  by  carrying  out  two 
experiments.  In  the  first  experiment  (phase  A) , substrate  is 
applied  at  a rate  in  excess  of  the  culture's  maximum  oxidizing 
rate.  By  measuring  SN  periodically,  the  term  dSN/dt  can  be 
found  and  used  to  calculate  UA  XB  A from  equation  (3-6)  . In  a 
subsequent  experiment  (phase  B) , substrate  is  applied  at  the 
appropriate  rate  so  that  Q2  SNo/V  is  approximately  equal  to 
one-half  UA  XB  A . The  resulting  steady-state  value  of  SN, 
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together  with  the  corresponding  UA  XB  A calculated  from 
equation  (3-7)  is  used  to  calculate  KNH  from  equation  (3-2) . 


Materials  and  Methods 


Summary  of  Experiments 

The  half-saturation  coefficients  for  ammonia  and  nitrite 
oxidization  were  determined  by  the  use  of  laboratory  scale 
fed-batch  experiments.  Two  2-L  glass  beakers,  two  stirring 
motors,  two  syringe  pumps  were  used  for  this  experiment.  A 
total  of  5 experiments  were  carried  out  for  each  oxidation 
step,  i.e.,  nitritif ication  and  nitratif ication. 

Microorganisms 

The  microorganisms  of  the  wasting  sludge  from  a bench- 
scale  contact  stabilization  system  were  used  in  this 
experiment.  In  order  to  maintain  same  SRT  in  contact 

stabilization  process,  a certain  amount  of  sludge  was 
withdrawn  from  reaeration  tank  through  a masterflex  pump  and 
stored  in  a refrigerator  to  keep  down  the  activity  of 
microorganisms.  This  refrigerated  sludge  was  were  aerated  to 
restore  the  activity  of  microorganisms  in  sludge. 

Experimental  Procedures 

This  experiment  consists  of  two  phases: 

i)  Phase  A was  to  estimate  the  unit  oxidization  rate 
with  excess  ammonia  or  nitrite  nitrogen. 
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ii)  Phase  B was  to  find  the  pseudo-steady-state 
concentration  of  either  ammonia  or  nitrite  nitrogen 
at  a selected  nitrification. 

The  value  of  UA  XB  A in  Equation  (3-6)  can  be  estimated 
based  on  the  substrate  feeding,  (Q  SN)/V,  and  accumulation 
rate,  dSN/dt,  in  Phase  A experiment.  UA  XB  A and  steady-state 
value  of  SN  can  be  also  evaluated  in  Phase  B experiment. 

800  mL  of  the  refrigerated  wasting  sludge  from  the 
contact  stabilization  system  was  aerated  for  at  least  one  hour 
to  acclimate  to  a room  temperature,  25±1°C,  before  inoculating 
two  2L  glass  beakers.  In  order  to  make  an  approximate 
concentration  of  270  mg/L  of  Mixed  Liquor  Suspended  Solids 
(MLSS ) in  each  reactor,  1300  mL  of  the  effluent  from  a 
secondary  clarifier  at  the  University  of  Florida  Sewage 
Treatment  Plant  (UFSTP)  was  added  to  each  beaker  with  400  mL 
of  acclimated  wasting  sludge.  Micronutrient  for  the  growth  of 
microorganisms  such  as  1.7  mL  of  28.7  g NaH2P04/L  , 1.7  mL  of 
0.06  g FeCl3- 6H20/L,  and  1.7  mL  of  22.5  g MgS04- 7H20/L  were 
added.  In  order  to  neutralize  the  acid  produced  from 
nitrification,  5 g of  sodium  bicarbonate  were  provided  to  each 
reactor. 

The  indoor  air  through  a pressure  reducer  (Dayton 
Electric,  Buffalo  Grove,  IL)  was  used  to  provide  the  air  to 
both  reactors.  The  debris  in  indoor  air  was  removed  through 
a cotton  filter  and  an  air  scrubber  was  used  to  humidify 
indoor  gas,  respectively.  A pressurized  carbon  dioxide  gas 
was  used  to  maintain  a pH  of  7 and  the  pure  oxygen  was 
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introduced  to  ensure  the  Dissolved  Oxygen  (DO)  concentration 
was  always  above  12  g 02/m3  during  experiment.  These  gases 
were  connected  to  indoor  air  tubes  through  Y connectors.  Two 
stone-to-air  diffusers  were  placed  at  the  bottom  of  each 
reactor  with  a Y-connector  as  in  Figure  3-1.  In  order  to 
maintain  a temperature  of  25±1  "C,  a water  bath  (Blue  Island, 
IL)  was  used.  Mixing  was  provided  using  a plastic  Nalgene 
propeller  attached  to  a stirring  motor. 

In  order  to  achieve  the  same  mixing  intensity  of  the 
contact  tank  in  contact  stabilization  system  in  each  batch 
reactor,  the  stirring  rate  of  the  mixer  in  each  batch  reactor 
was  adjusted  using  the  following  relationship  (American  Water 
Works  Association  Research  Foundation,  1991) : 


N, 


3/2 


nb3/2 


V, 


1/2 


V, 


1/2 


(3-8) 


'C  ’B 

where,  Nc  = revolution  rate  of  the  mixer  in  contact  tank,  rpm 

NB  = revolution  rate  of  the  mixer  in  each  batch 
reactor,  rpm 

Vc  = volume  of  contact  tank,  L 

VB  = volume  of  each  batch  reactor,  L 


From  equation  (3-8) , Ns  is 


Vb/Z 

VcZ 


N*'2 


2/3 


(3-9) 
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Ng  is  to  be  685  rpm  since  Vg,  Vc  and  Nc  are  1.7  L,  2.5  L and 
780  rpm,  respectively.  This  revolution  rate  was  maintained  in 
both  batch  reactors  during  experiment.  Digital  Optical 
Tachometer  (Model-IV,  Monarch  Instrument,  Amherst,  NH)  was 
used  to  measure  the  rpm  of  the  mixer  in  each  tank. 

In  the  experiment  of  the  half-saturation  coefficient  for 
Nitrosomonas . a 5000  g NH4+1-N/m3  ammonium  chloride  solution 
was  fed  into  the  Phase  A reactor  by  a syringe  pump  (model  355, 
Orion,  Boston,  MA)  through  a length  of  Tygon  tube  (0.78  mm 
I.D.)  which  reached  to  the  bottom  of  a reactor.  The  feed  rate 
to  phase  A reactor  was  3.47  mL/hr  to  make  a rate  in  excess  of 
the  culture's  maximum  nitrification  rate.  This  rate  was 
selected  by  trial  and  error.  Initially,  5 mL  of  the  above 
solution  was  spiked  to  make  a 15.0  g NH4+1-N/m3  in  the  Phase  A 
reactor. 

In  the  Phase  B experiment,  the  procedure  of  Phase  A 
experiment  was  used  except  for  the  feed  rate  of  a ammonium 
chloride  solution  and  the  initial  spike.  The  feed  rate,  0.68 
mL/hr,  was  selected  to  make  a pseudo-steady-state  condition  by 
trial  and  error.  The  reactor  was  spiked  to  make  5 g NH4+1-N/m3 
by  1.7  mL  of  a 5000  g NH4+1-N/m3  ammonium  chloride  solution. 

In  experiments  for  the  half-saturation  coefficient  for 
Nitrobacter . a 5000  g NO^-N/m3  sodium  nitrite  solution  was 
fed  into  the  phase  A reactor.  Since  a high  concentration  of 
nitrite,  i.e.,  30  g NO^-N/m3  showed  the  substrate  inhibition 
to  Nitrobacter  as  in  Figure  3-2,  no  nitrite  was  added  to 
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phase  A reactor  for  initial  2 hours  experimental  time.  A 
feeding  concentration  in  Figure  3-2  represents  the  estimated 
N02'-N  concentration  based  on  feeding  rate  without  oxidization 
reaction.  During  this  period,  microorganisms  were  accumulated 
to  nitrite.  2 mL  of  a 5000  g NO^-N/m3  sodium  nitrite  solution 
initially  was  spiked  to  make  6 g N02-N/m3  in  each  reactor. 

In  Figure  3-2,  a feeding  concentration  refers  the  estimated 
concentration  of  nitrite  nitrogen  based  upon  the  sodium 
nitrite  feeding  rate. 

In  the  phase  B experiment,  a feed  rate  of  0.63  mL/hr  of 
5000  g NO^-N/m3  sodium  nitrite  solution  was  chosen  to  make  a 
pseudo-steady-state  by  trial  and  error.  1 mL  of  5000  g N02-1- 
N/m3  sodium  nitrite  solution  was  initially  spiked  to  bring  the 
concentration  to  3 g N02'1-N/m3  in  the  phase  B reactor. 

Analytical  Methods 

Ammonia  Nitrogen.  Ammonia  nitrogen  was  analyzed  using 
a gas  sensing  combination  electrode  (Orion  model  9512, 
Boston,  MA)  with  a Digital  Ionanalyzer  (Orion  model  701A, 
boston,  MA)  according  to  417  E in  the  Standard  Methods  (APHA, 
et  al.,  1985).  30  mL  of  samples  from  each  reactor  taken  by  a 
syringe  was  filtered  through  a pore  size  0.45  yum  membrane 
filter  (Gelman  GN-6,  Ann  Arbor,  MI)  after  prefiltration  with 
a pore  size  1.6  /xm  glass  microfibre  filter  (Whatman  GF/C, 
Maid  stone,  England) . 


N02-N  CONC.,  mg/L 


29 


Measured  Cone.  — ■ — Feeding  Cone. 


Figure  3-2.  Nitrite  nitrogen  inhibition  to  Nitrobacter . 
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5 ml  of  filtered  samples  were  used  to  analyze  the 
concentrations  of  ammonia  nitrogen.  In  order  to  minimize  the 
ammonia  volatilization,  a 17  ml  glass  vial  having  a very 
narrow  neck  and  therefore  a very  small  surface  exposed  to  the 
air  were  used.  The  mixing  was  provided  by  a magnetic  stirrer 
(Model  PC353 , Corning,  Corning,  NY)  with  a micro  bar  (01/8"  x 
Ll/2")  before  adding  0.3  ml  of  a 5M  NaOH  solution.  This  NaOH 
solution  contained  0.05M  Disodium  Ethylene  Diamine  Tetra 
Acetic  Acid  (EDTA) , 10%  (V/V)  methanol  and  5M  NaOH.  After  a 
stable  reading  was  displayed,  the  mV  value  was  recorded. 
After  calibration  against  standard  solutions,  the 
concentrations  of  ammonia  nitrogen  in  samples  were  determined. 

Mixed  Liquor  Suspended  Solids  fMLSS) . Samples  of  5 mL 
from  each  reactor  and  100  mL  from  an  effluent  tank  were  taken 
to  measure  the  concentrations  of  Mixed  Liquor  Suspended  Solids 
(MLSS ) . These  samples  were  filtered  through  a predried  and 
preweighted  a 1.2  fim  glass  microfibre  filter  (Whatman  GF/C, 
Maidstone,  England) . After  drying  in  an  oven  at  a temperature 
of  103 "C  for  at  least  three  hours,  the  MLSS  were  estimated  by 
the  difference  of  filtration  paper  weight. 

Nitrite  and  Nitrate  Nitrogens.  Two  oxidized  nitrogen 
forms,  nitrite  and  nitrate  nitrogen  were  measured  using  an  Ion 
Chromatography  (Dionex  model  2 000i,  Sunnyvale,  CA)  . 
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The  extra  filtered  samples,  approximately  5 mL,  after  an 
ammonia  nitrogen  analysis  were  used  to  inject  into  this 
instrument.  Since  a sample  had  high  concentrations  of  nitrate 
nitrogen,  samples  were  diluted  to  a factor  of  20  times.  A 
standard  solution  having  11.288  mg  N02‘-N/L  and  15.215  mg  N03'- 
N/L  was  used  to  calibrate. 


Results  of  Experiments 

Experimental  results  of  the  half-saturation  coefficients 
for  ammonia  and  nitrite  oxidizers  are  shown  in  Tables  3-1  and 
3-2,  respectively.  In  these  Tables,  A and  B refers  Phase  A 
and  B experiments,  respectively.  Phase  A was  used  to  obtain 
uAXB  A and  UAXB  A and  a steady-state  value  were  estimated  at 
Phase  B.  The  half-saturation  coefficient  for  Nitrosomonas 
ranged  from  0.728  g NH4+1-N/m3  to  0.455  g NH4+1-N/m3  with  an 
average  of  0.58  g NH4+1-N/m3  at  a temperature  of  25°C  and  an 
average  pH  of  7.  The  duration  time  for  this  experiment  was  7 
hours.  The  results  of  5 pairs  of  experiments  are  shown  in 
Figures  3-3,  3-4  and  3-5. 

The  values  for  the  half-saturation  coefficient  for 
Nitrobacter  ranged  from  0.883  g N02'1-N/m3  to  0.516  g NO^-N 
with  an  average  of  0.719  g NO^-N/m3  under  same  conditions  of 
the  experiment  for  Nitrosomonas . The  results  of  5 pairs  of 
experiments  are  shown  in  Figures  3-6,  3-7  and  3-8. 


Table  3-1.  Summary  of  KNH4A  calculation  for  the  experimental  data. 
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Table  3-2.  Summary  of  KN02A  calculation  for  the  experimental  data. 
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Figure  3-3.  Determination  of  KNH4A  using  fed  batch  experiments 
7/9  A,  7/9  B,  7/11  A and  7/11  B (■:  Data  used 
for  estimation)  . 
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Figure  3-4.  Determination  of  KNHAA  using  fed  batch  experiments 
7/12  A,  7/12  B,  7/17  A and  7/17  B (■:  Data  used 
for  estimation) . 
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Figure  3-5.  Determination  of  KNH,.  using  fed  batch  experiments 
7/18  A and  7/18  B (■:  Data  used  for  estimation). 
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Figure  3-6.  Determination  of  KN02A  using  fed  batch  experiments 
10/25  A,  10/25  B,  10/28  A and  10/28  B (■:  Data 
used  for  estimation) . 
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Figure  3-7.  Determination  of  KN02A  using  fed  batch  experiments 
10/31  A,  10/31  B,  11/1  A and  11/1  B (■:  Data 
used  for  estimation) . 
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Determination  of  KN02.  using  fed  batch  experiments 
11/8  A and  11/8  B (■:  Data  used  for  estimation)  . 
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CHAPTER  4 


DETERMINATION  OF  THE  EFFECTIVE  MAXIMUM  SPECIFIC 
GROWTH  RATES  OF  NITROGEN  OXIDIZING  BACTERIA1 

Introduction 

The  maximum  specific  growth  rates  of  the  nitrogen 
oxidizing  bacteria  in  activated  sludge  are  most  crucial 
parameters  describing  the  growth  of  the  autotrophic  bacteria. 
Two  maximum  specific  growth  rates  are  Aanh4  and  Aano2  for  ammonia 
oxidizing  and  nitrite  oxidizing  bacteria,  respectively.  In 
this  chapter,  Aa  refers  to  either  microorganism's  maximum 
specific  growth  rates.  The  correct  estimation  of  Aa  will 
determine  the  nitrification  capacities  in  a wastewater 
treatment  plant,  the  possibility  of  nitrite  build-up  in  an 
aeration  tank  and  the  minimum  Sludge  Retention  Time  (SRT) , 
below  which  washout  of  the  nitrogen  oxidizers  would  occur. 

Since  the  Aa  has  a big  effect  on  nitrification,  the 
estimation  of  these  parameters  is  very  important  for 
simulating  nitrification  using  a computer  model. 


1 Nitrogen  oxidizing  bacteria  include  ammonia  and  nitrite 
oxidizing  bacteria 
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Basis  for  the  Determination  of  flA  - bA  Using  a Batch  Culture 
When  nitrogen  oxidizing  bacteria  are  grown  in  a batch 
reactor  at  non-limiting  nitrogen  and  dissolved  oxygen 
concentrations,  their  concentrations  vary  according  to  the 
following  equation: 

dXB  A 

 *  1 * = (Aa  - bA  ) X (4-1) 

dt 

Xg  A = active  autotrophic  biomass,  g/m3 

= XBANH4  in  ammonia  oxidization 

= XBAN02  in  nitrite  oxidization 

bA  = autotrophic  decay  coefficient  in  either  ammonia 
or  nitrite  bacteria 

— = growth  rate  of  autotrophic  biomass,  g/m3 

dt 


where , 


^B,A 


If  the  production  of  nitrite  is  completely  growth  associated 
and  nitrite  were  not  further  oxidized  to  nitrate  in  the  case 
of  ammonia  nitrogen,  we  would  have: 


d[N02'1-N] 

dt 


1 

Aanh4  XBANH4 
Y 

■‘•ANHA 


(4-2) 


where, 


dCNO^-N] 

dt 


accumulation  rate  of  nitrite, 
g/m3  day 


Y 


ANH4 


= growth  yield  for  ammonia  oxidizers,  g/g 
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Equation  (4-1)  can  be  integrated  subject  to  the  initial 
condition  XBANH,(0)  = XBANH4o  to  give: 

(AANH4  ~ t 

XBANH4  = XBANH4o  e (4“3) 

combining  equations  (4-2)  and  (4-3) , the  following  expression 
is  obtained: 


d[N02'1-N]  1 


Aaiju/  X, 


(AaNH4  " bA>  t 


dt 


ANH4  BANH4o 


‘ ANH4 


(4-4) 


which  upon  integration  yields: 


[N02"  -N]t  = 


AaNH4  XBANH4o  e 
YANH4  (AANH4  “ bA^ 


(AanH4  b 


+ [N02’1-N] Q - 


AaNH4  XBANH4o 


YANH4  ^ AANH4  bA) 


(4-5) 


where  [N02'1-N]0  is  the  initial  nitrite  nitrogen  concentration. 
With  the  exception  of  an  initial  start-up  time,  the  sum  of  the 
second  and  third  terms  in  the  right  hand  side  of  equation  (4- 
5)  will  be  negligible  since  they  are  constants  while  the  first 
term  increases  exponentially. 

If  the  natural  logarithm  is  taken  in  both  sides  without 
the  second  and  third  terms  in  equation  (4-5) , the  following 
equation  can  be  derived: 

, AAnh4  XBANH4o 

In  [N02"  -N] t = In  + (AANH4-bA)*t  (4-6) 

YANH4(AAnh4  “ bA  ^ 
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In  the  mixed  culture  experiments  carried  out,  most  of 
nitrite  produced  was  further  oxidized  to  nitrate  by  nitrifying 
bacteria.  The  sum  of  the  nitrite  and  nitrate  nitrogen 

concentrations  thus  gives  the  total  amount  of  nitrite  that 
would  be  present  in  the  absence  of  the  nitrite  oxidizers. 
Assimilation  of  nitrite  or  nitrate  as  a nitrogen  source  for 
biosynthesis  can  be  neglected,  since  ammonia  was  present  in 
excess  during  experimental  period  and  is  the  preferred 
nitrogen  source  for  microorganisms.  Consequently,  plotting 
ln([N02'1-N]  + [NOj^-N])  versus  time  gives  a slope  of  (Aanh4  ~ 
bA)  which  is  the  effective  specific  growth  rate  of 
Nitrosomonas . Due  to  the  terms  neglected  in  equation  (4-6) , 
a deviation  from  straight  line  in  the  plot  is  expected  for  the 
first  few  experimental  data. 

For  nitrite  oxidation,  the  following  equation  can  be 
derived: 

. ^AN02  ^BAN02o 

In [N03‘  -N] t = In  + (£AN02-bA)  t (4-7) 

^AN02  (^AN02  ” bA^ 

where,  [NO^-N^  = concentration  of  nitrite  nitrogen  at  t, 

g/m3 

Yano2  = growth  yield  for  nitrite  oxidizers,  g/g 

By  plotting  ln[N03'1-N]  versus  time,  a slope  (AAN02-bA) 
which  is  the  effective  specific  growth  rate  of  Nitrobacter 


could  be  estimated. 
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Materials  and  Methods 


Summary  of  Experiments 

The  specific  maximum  growth  rates  for  Nitrosomonas  and 
Nitrobacter  were  determined  using  laboratory  fed-batch 
reactors  as  in  the  experiments  for  the  half-saturation 
coefficients.  A total  of  six  experiments  each  for 

Nitrosomonas  and  for  Nitrobacter  were  carried  out. 

Microorganisms 

The  microorganisms  described  in  Chapter  3 were  used  in 
these  experiments. 

Experimental  Procedures 

A refrigerated  wasting  sludge  was  aerated  for  at  least  1 
hour  to  acclimate  to  room  temperature  of  25 °C.  170  ml  of  the 

acclimated  wasting  sludge  were  added  to  two  2 L beakers  in  a 
water  bath  for  temperature  control  along  with  micronutrients 
and  1.53  L of  secondary  clarifier  effluent  from  the  UFSTP.  In 
order  to  neutralize  acids  produced  from  nitrification,  5 g of 
sodium  bicarbonate  were  added  to  both  beakers.  The 
experimental  apparatus  was  the  same  as  in  Figure  3-1.  In 
order  to  avoid  lack  of  substrate,  i.e.,  ammonia  nitrogen,  and 
inhibition  of  nitrification  due  to  high  concentrations  of  Free 
Ammonia  (FA) , a concentration  of  approximately  10  g NH4+-N/m3 
was  maintained  during  experiments.  At  the  beginning  of  each 
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experiment,  ammonia  nitrogen  was  provided  by  spikes.  After  the 
estimation  of  a conversion  rate,  a syringe  pump  was  used.  A 
solution  of  ammonium  sulfate,  5000  g NH4+-N/m3,  was  used  as  a 
source  of  ammonia  nitrogen  instead  of  ammonium  chloride  since 
the  peak  of  nitrite  overlapped  with  that  of  chloride  in  the 
printout  of  the  ion  chromatograph,  whereas  a 5000  g N02‘-N/m3 
sodium  nitrite  solution  was  used  as  the  source  for  nitrite 
oxidizing  bacteria.  Micronutrients  were  added  as  described  in 
chapter  3.  At  the  end  of  each  experiment,  the  concentrations 
of  MLSS  in  both  reactors  were  measured. 

Analytical  Methods 

Ammonia  Nitrogen 

The  method  described  in  Chapter  3 was  used. 

Nitrate  and  nitrite  nitrogen 

The  method  described  in  Chapter  3 was  used. 

Mixed  Liquor  Suspended  Solids 

The  method  described  in  Chapter  3 was  used. 

Results  of  Experiments 

Experimental  results  for  ammonia  and  nitrite  nitrogen 
oxidation  are  shown  in  Tables  4-1  and  4-2,  respectively.  The 
values  for  the  effective  maximum  specific  growth  rate  of 
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Nitrosomonas  ranged  from  0.6799  day'1  to  0.5819  day'1  with  an 
average  of  0.6510  day'1  under  conditions  of  a temperature  of 
25 °C  and  an  average  pH  of  7.  The  duration  time  for  the 
experimental  phase  was  about  3 days  while  the  production  of 
oxidized  nitrogen  was  360  gN/m3.  The  results  of  six 
experiments  are  shown  in  Figures  4-1  and  4-2. 

Values  for  the  effective  maximum  specific  growth  rate  of 
Nitrobacter  ranged  from  0.8741  day'1  to  0.7  623  day'1  with  an 
average  of  0.8204  day'1  under  the  similar  same  conditions  as 
for  the  Nitrosomonas  experiments.  The  exponential  phase 
lasted  for  about  3 days  while  production  of  nitrate  nitrogen 
was  about  3 00  gN/m3.  The  results  of  six  experiments  are  shown 
in  Figures  4-3  and  4-4. 

As  seen  from  semi  logarithmic  plots  (Figures  4-1  - 4)  for 
all  the  experiments,  the  majority  of  data  fit  straight  lines, 
as  predicted  from  the  theory.  In  determining  (£A  - bA)  , only 
the  data  from  the  exponential  phase,  i.e.,  those  lying  on  a 
straight  line  were  used.  For  all  experiments,  the  data  which 
were  used  to  compute  the  effective  maximum  specific  growth 
rates  gave  very  good  correlation  coefficients,  r2  > 0.99. 
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Table  4-1.  Summary  of  experimental  results  for  the  effective 
specific  growth  rate  of  Nitrosomonas . 


Exp.  # 

#ANH4‘bA 

1/d 

r2 

Data 

Points 

pH 

Avg.  Temp.”  C 

MLSS 

9/m3 

1 

0.6705 

0.992 

10 

6.8-72 

25 

200.5 

2 

0.6482 

0.992 

7 

6.8-7.2 

25 

200.8 

3 

0.5819 

0.993 

8 

6.8-7.2 

25 

195.6 

4 

0.6799 

0.998 

9 

6. 8-7,2 

25 

195.6 

5 

0.6663 

0.995 

10 

68-7.2 

25 

210.5 

6 

0.6593 

0.996 

10 

6. 8-7. 2 

25 

210.6 

Average  Value  of  AANH4  ~ t>A  = °-6510  i/d 
Standard  Deviation  = 0.0408  1/d 
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Table  4-2.  Summary  of  experimental  results  for  the  effective 
specific  growth  rate  of  Nitrobacter . 


Exp.  # 

^AN02"^A 

1/d 

r2 

Data 

Points 

pH 

Avg.  Temp.°C 

MLSS 

9/m3 

1 

0.8169 

0.998 

9 

6. 8-7.2 

25 

189.4 

2 

0.8715 

0.999 

9 

6. 8-7.2 

25 

189.4 

3 

0.8088 

0.994 

8 

6. 8-7.2 

25 

185.5 

4 

0.8741 

0.997 

8 

6.8-7.2 

25 

185.6 

5 

0.7623 

0.991 

8 

6. 8-7.2 

25 

187.8 

6 

0.7890 

0.998 

9 

6.8-7.2 

25 

188.3 

Average  Value  of  Aan02  “ = °*8204  1/d 

Standard  Deviation  = 0.0324  1/d 


In  (Oxidized  N Cone  ),  mg/L  Ln(Oxidized  N Cone),  mg/L 
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Figure  4-1.  Nitritif ication  results  from  experiments  1,2,3  and 
4 (■:  Data  used  for  estimation) . 
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Figure  4-2.  Nitritif ication  results  from  experiments  5 and  6 
(■:  Data  used  for  estimation) . 
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Figure  4-3 


Nitratif ication  results  from  experiments  1,2,3 
and  4 (■:  Data  used  for  estimation) . 


Ln(N03-N  Cone),  mg/L  Ln(NQ3^  Cone),  mg/L 
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Figure  4- 
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Avg.  pH  = 7.0 
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4.  Nitratif ication  results  from  experiments  5 and  6 
(I:  Data  used  for  estimation) . 


CHAPTER  5 


CONTROL  OF  NITRIFICATION 
IN  CONTACT  STABILIZATION  PROCESS 

Introduction 

The  contact  stabilization  process  is  a modified 
activated  sludge  system  that  has  been  used  to  treat  domestic 
wastewater  with  less  aeration  tank  volume  compared  to  common 
activated  sludge  systems  (Metcalf  & Eddy,  1991)  . Although 
partial  nitrification  in  the  contact  stabilization  process 
was  studied  by  several  researchers  (Gujer  and  Jenkins,  1975; 
Alexander,  et  al.,  1980;  Dold,  et  al.,  1980),  control 
strategies  for  nitrification  were  not  fully  investigated. 

In  this  study,  partial  nitrification  was  evaluated  by 
changing  two  operation  parameters;  Sludge  Retention  Time 
(SRT)  and  Sludge  Recycle  Ratio.  At  SRTs  of  4.8  days  and  6 
days,  three  different  sludge  recycle  ratios,  i.e.,  50%,  80% 
and  110%  of  influent  flow  rate  (0.0125  m3/day)  , were  used  to 
find  their  effects  on  nitrification.  In  order  to  investigate 
the  effects  of  sludge  recycle  ratio  further  on  nitrification, 
20%  and  140%  of  recycle  ratios  were  added  at  6 days  of  SRT. 
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Materials  and  Methods 


Microorganisms 

A mixed  culture  taken  from  a contact  tank  at  the 
University  of  Florida  Sewage  Treatment  Plant  (UFSTP) , 
Gainesville,  Florida,  was  used.  A phytoplankton  screen,  mesh 
size  3 65  with  an  opening  diameter  of  20  cm,  was  used  to 
remove  sand  before  inoculation,  since  a small  diameter  (1/4") 
Tygon  tube  was  used  in  a bench-scale  contact  stabilization 
system. 

Synthetic  Wastewater 

The  wastewater  in  this  study  was  assumed  to  have  two 
waste  streams:  the  effluent  from  a primary  activated  sludge 
system  and  spent  nutrient  solution  from  a hydroponic  culture. 
Since  a primary  activated  sludge  is  assumed  to  carry  out  only 
the  oxidation  of  organic  materials  without  nitrification, 
the  effluent  has  high  concentrations  of  ammonia  nitrogen. 
Awong  (1990)  reported  the  spent  nutrient  solution  had  several 
sugar  components,  i.e.,  Glucose,  Fructose,  and  Sucrose 
(Department  of  Environmental  Engineering  Sciences,  et  al, 
1990) . 

In  order  to  produce  the  same  constituents  in  the 
synthetic  wastewater,  ammonia  nitrogen  and  sugar  components 
were  added  to  the  effluent  from  the  University  of  Florida 
Sewage  Treatment  Plant  (UFSTP) . The  synthetic  wastewater 
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consisted  of  Chemical  Oxygen  Demand  (COD)  constituents, 
ammonia  nitrogen  components  and  three  micronutrients.  Table 
5-1  shows  the  concentration  of  stock  solutions  and  the  added 
amounts  to  0.075  m3  of  effluent  to  make  the  target  wastewater 
composition.  Table  5-2  shows  the  constituent  concentrations 
in  synthetic  wastewater  after  adding  all  components. 


Experimental  Procedures 

A Contact  Stabilization  System  was  used  to  investigate 
nitrification  in  biological  processes.  This  system  had  two 
bioreactors,  the  contact  and  reaeration  tanks  as  shown  in 
Figure  5-1.  Two  DO  (Dissolved  Oxygen)  controller  complexes 
were  used  to  maintain  DO  concentrations  at  4.0±0.5  g02/m3  in 
both  reactors.  DO  electrode/probe  (900  series,  New 
Brunswick,  Edison,  NJ) , DO  analyzer  (Model  DO-40  New 
Brunswick,  Edison,  NJ) , and  recorder/controller  (Model  305- 
720,  Gulton  Rustrack,  Manchester,  NH)  were  included  in  a DO 
controller  complex. 

The  influent  and  a sludge  recycle  pump  had  computerized 
drive  motors  (Model  7550-90,  Cole  Parmer,  Chicago,  IL)  and 
they  were  controlled  by  an  IBM-compatible  computer. 
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Table  5-1.  Constituents  of  stock  solutions. 


Constituent 

Concentration 
in  stock  solution 

Amount 
per  75L  in  SE 

Fructose 

81.3  g/L 

50.0  mL 

Glucose 

81.3  g/L 

50.0  mL 

Sugar 

77.6  g/L 

50.0  mL 

Sodium  Acetate 

277.0  g/L 

lOO.OmL 

Ammonium  Chloride 

319.0  g/L 

lOO.OmL 

Ammonium  Sulfate 

395.0  g/L 

lOO.OmL 

Magnesium  Sulfate 

22.5  g/L 

75.0  mL 

Ferric  Chloride 

0.25  g/L 

75.0  mL 

Sodium  Phosphate 

33.0  g/L 

75.0  mL 

Sodium  Carbonate 
(Baking  Soda) 

20.0  g 

NOTE:  SE  = THE  SECONDARY  CLARIFIER  EFFLUENT  AT  THE  UNIVERSITY  OF  FLORIDA 


SEWAGE  TREATMENT  PLANT 
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Table  5-2 . 


Concentrations  of  constituents  in  synthetic 
wastewater. 


Constituents 

Concentration 

Unit 

Biodegradable 
Soluble  COD 

347.9 

mg  COD/ 

Unbiodegradable 
Soluble  COD 

38.5 

mg  COD/ 

Ammonia  Nitrogen 

221.5 

mg  N/L 

Nitrate  Nitrogen 

12.2 

mg  N/L 

Alkalinity  as  CaC03 

2600.0 

mg/L 

REFRIGERATOR 
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Each  motor  had  a Masterflex  L/S  standard  pump  head  (Model 
7016,  Cole  Parmer,  Chicago,  IL)  . A Masterflex  pump  (Model 
7553-10,  Cole  Parmer,  Chicago,  IL)  with  a head  (Model  7016, 
Cole  Parmer,  Chicago,  IL)  was  used  for  a sludge  wasting  pump 
to  maintain  a constant  SRT. 

A Norprene  tube  was  used  for  the  influent  and  sludge 
recycle  pumps  since  this  tube  could  last  for  more  than  5000 
hour  operation  time.  The  sludge  wasting  pump  used  a silicone 
tube  since  this  pump  was  not  operated  continuously.  The 
sludge  wasting  pump  was  operated  for  2 minutes  each  hour 
using  a timer. 

A solenoid  valve  (Grainger  Model  741150115,  Arlington 
Height,  IL)  was  used  to  control  the  flow  of  oxygen  from  an 
oxygen  cylinder  through  a rotameter  (Model  N062-01,  Cole 
Parmer,  Chicago,  IL)  since  the  Oxygen  Uptake  Rate  (OUR)  was 
too  high  in  the  contact  tank  to  maintain  a target 
concentration  of  DO  by  aeration  only.  When  a fine  bubble 
diffuser  (Kondon  mist-air  62503,  Hayward,  CA)  was  placed  at 
the  bottom  of  reactor,  there  was  interference  of  the  DO  probe 
due  to  attachment  of  oxygen  bubbles  to  the  probe  membrane. 
In  order  to  solve  this  problem,  a fine  bubble  stone  diffuser 
along  with  a coarse  bubble  diffuser  (Kondon  mist-air  62501, 
Hayward,  CA)  was  placed  at  1/4  depth  from  the  bottom.  A 
coarse  air  bubble  diffuser  provided  extra  mixing  to  the 
contact  tank.  Air  was  supplied  continuously  by  an  air  pump 
(Willinger  Bros.,  Inc.,  Fort  Lee,  NJ) . 
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In  order  to  maintain  a pH  of  7.0  ± 0.2,  carbon  dioxide 
was  supplied  to  the  contact  tank  through  a rotameter  (Model 
N042-15,  Cole  Parmer,  Chicago,  IL) . The  carbon  dioxide  tube 
was  inserted  into  the  air  supply  tube  through  a needle. 

A water  bath  and  circulator  ( Lauda/ Brinkmann,  Westtowry, 
NY)  was  used  to  maintain  a temperature  of  25  ± 1 °C  in  both 
reactors.  The  water  was  circulated  through  1/4"  glass  tubes 
located  along  the  wall  in  each  reactor. 

The  contact  tank  had  a mixer  (Model  65748,  Precision 
scientific,  Chicago,  IL)  with  a Nalgene  polyethylene  stirrer 
operated  at  780  rpm.  An  impeller  was  located  at  the  bottom 
of  the  reactor.  A mixer  (Dayton  Electric  manufacture  Co., 
Chicago,  IL)  with  two  pairs  of  flat  type  impellers  rotating 
at  460  rpm  was  used  to  mix  the  contents  of  the  reaeration 
tank.  One  pair  was  located  at  the  middle  and  the  other  was 
placed  at  the  bottom  of  the  reactor,  respectively. 

The  sludge  transfer  tube,  1/2"  of  diameter,  between  the 
reaeration  tank  and  the  contact  tank  was  aerated  through  a 
needle  to  avoid  anoxic  conditions.  Air  was  provided  by  an 
air  pump  (Willinger  Bros.,  Inc.,  Fort  Lee,  NJ)  . Since  a 
computerized  system  required  a reliable  electric  power 
source,  a battery  backup  system  (Model  BC-750,  Tripplite, 
Chicago,  IL)  was  used. 

The  influent  pump  with  a flow  rate  0.025m3/day  was  only 
operated  for  30  minutes  in  every  hour  since  a continuous 
inflow  was  found  to  cause  sludge  bulking  in  the  contact  tank. 
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The  sludge  recycle  pump  was  operated  with  different  flow 
rates  depending  on  the  target  recycle  ratio.  At  a ratio  of 
110  percent,  the  flow  rate  was  9.54  ml/min.  Flow  rates  of 
6.94  ml/min  and  4.34  ml/min  were  used  for  ratios  of  80  and  50 
percent,  respectively.  At  recycle  rates  of  20  and  140 
percents  flow  rates  were  1.73  and  12.14  mL/min,  respectively. 

The  flow  rate  of  the  sludge  wasting  pump  was  adjusted 
according  to  the  target  SRT,  after  measurement  of  the 
concentrations  of  Mixed  Liquor  Suspended  Solids  (MLSS)  in 
each  reactor  and  the  effluent. 


Analytical  Methods 

Samples  from  each  reactor  were  taken  to  analyze  the 
concentrations  of  ammonia,  nitrite  and  nitrate  nitrogen 
daily.  Concentrations  of  MLSS  were  also  measured  in  each 
reactor  and  the  effluent  daily.  Soluble  Chemical  Oxygen 
Demand  (SCOD)  and  Soluble  Total  Kjeldahl  Nitrogen  (STKN)  in 
both  reactors  through  a 0.45  jum  membrane  filter  was  measured 
weekly. 

Ammonia  Nitrogen  Ammonia  nitrogen  was  analyzed  using 
the  Method  described  in  Chapter  2 . 
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Nitrite  and  Nitrate  Nitrogen  Two  oxidized  nitrogen 
forms,  nitrite  and  nitrate  nitrogen  were  measured  using  Ion 
Chromatograph  (Dionex  model  2000i,  Sunnyvale,  CA)  described 
as  in  Chapter  3 . 

Mixed  Liquor  Suspended  Solids  (MLSS)  The  method 
described  in  Chapter  3 was  used  to  measure  the  concentration 
of  Mixed  Liquor  Suspended  Solids  (MLSS) . 

Chemical  Oxygen  Demand  Chemical  Oxygen  Demand  (COD) 
was  analyzed  using  a COD  reactor  and  COD  vials  (Model  45600C, 
HACH,  Loveland,  CO)  . A 2 mL  sample  filtered,  as  in  the 
ammonia  nitrogen  analysis,  was  introduced  to  a vial.  After  2 
h digestion  in  a COD  reactor,  a digested  sample  was  titrated 
against  a 0.0125  M ferrous  ammonium  sulfate  (FAS)  solution. 
This  FAS  solution  was  standardized  against  a 0.025  M 
potassium  dichromate  standard  solution  before  the  titration. 
In  case  of  total  COD  estimate  of  a mixed  culture,  a blender 
was  used  to  homogenize  the  content  before  the  digestion.  In 
order  to  avoid  the  COD  resulting  from  nitrite  in  a sample,  an 
appropriate  amount  of  sulfamic  acid  was  added  to  a COD  vial 
when  a sample  had  nitrite  (American  Public  Health  Association 
et  al,  1985)  . 

Total  Kieldahl  Nitrogen  Total  Kjeldahl  Nitrogen 
(TKN)  was  measured  according  to  NO.  00625  in  Methods  for 
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Chemical  Analysis  of  Water  and  Wastes  (U.S.  EPA,  1979). 

The  samples  for  soluble  TKN  were  filtered  as  in  the  ammonia 
nitrogen  analysis  except  for  a filterate  nitrate  removal 
procedure.  Since  higher  nitrate  concentrations  (10  times  or 
higher  than  the  TKN  level)  can  result  in  low  TKN  values,  an 
anion  exchange  resin  (Rohm  and  Hass  Co.  , Philadelphia,  PA) 
was  used  to  remove  nitrates  in  the  filtrate.  The  reguired 
amount  of  resin  was  calculated  according  to  Hutton  (1991) . 

In  order  to  avoid  nitrate  interference  in  total  TKN 
measurement,  separation  of  solids  by  a centrifuge  was  adopted 
to  measure  TKN  of  solids  since  an  anion  exchange  resin  cannot 
be  separated  from  solids  after  nitrate  removal  from  samples. 
After  30  minutes  of  centrifugation,  the  supernatant  having 
nitrate  was  discarded  using  a Pasteur  pipet  and  ammonia-free 
water  was  filled  to  make  the  original  volume  (20  mL)  . A 
block  digester  (Model  BD-40,  Technicon,  Tarr  Town,  NY)  with 
75  mL  constricted  digester  tubes  (Bran+Luebbe,  Elmsford,  NY) 
was  used  to  digest  samples  along  with  a digest  solution 
(Sulfuric  acid-mercuric  sulfate-potassium  sulfate  solution) 
for  1 hour  at  160 °C  and  2 and  half  hours  at  360 °C.  Organic 
nitrogen  in  sample  was  to  convert  to  ammonia  nitrogen  after 
digestion  procedure.  Ammonia  from  organic  nitrogen  in  a 
digested  solution  and  inorganic  ammonia  nitrogen  was  measured 
using  same  method  described  in  ammonia  nitrogen  measurement. 
At  this  analysis,  a NaOH  solution,  however,  had  0.3  mL  of  a 
10  N NaOH  solution  including  Nal  and  EDTA. 
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Results  of  Experiments 

The  steady-state  condition  was  determined  after 
reviewing  a graph  having  six  components,  i.e.,  ammonia 
nitrogen,  nitrite  nitrogen,  nitrate  nitrogen  and  the 
concentrations  of  MLSS  in  each  reactor  at  different  date. 
Figure  5-2  shows  the  experimental  data  under  operating 
conditions  at  4.8  days  of  SRT  with  80%  of  sludge  recycle 
ratio. 

The  steady-state  values  for  this  experiment  were 
summarized  in  Table  5-3.  The  difficulty  in  maintaining  a 
constant  temperature  (25°C)  in  both  reactors  during  operation 
caused  high  concentrations  of  suspended  solids  in  the 
effluent  and  occasionally  lengthened  the  time  to  reach  a 
steady-state  under  given  conditions. 


~i/m  Rin  ' nnn^  nnRmiiK 
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Figure  5-2.  Experiment  data  at  4.8  days  of  SRT  with  80%  of 
sludge  recycle  ratio 


MLSS  Cone.,  mg/L 
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CHAPTER  6 

MODEL  DEVELOPMENT  AND  COMPUTER  SIMULATION 
Introduction 

The  International  Association  on  Water  Pollution  Research 
and  Control  (IAWPRC)  task  group  on  biological  modeling 
developed  a model  for  single  sludge  systems  that  accounted  for 
eight  process  rates  and  thirteen  wastewater  and  biological 
components  (Henze  et  al , 1986).  Their  model  divides  active 
biomass  into  two  components  (heterotrophic  and  autotrophic) 
and  nitrogen  into  four  components  (ammonia  nitrogen, 
particulate  degradable  organic  nitrogen,  soluble  degradable 
organic  nitrogen,  nitrate  nitrogen) . Two  additional  nitrogen 
species,  particulate  inert  organic  nitrogen  and  soluble  inert 
nitrogen  are  recognized  in  the  model  but  they  do  not  enter  the 
reaction  equations. 

In  this  study,  the  heterotrophic  and  autotrophic  biomass 
components  were  divided  into  two  groups,  respectively,  so  that 
the  model  contains  four  biomass  components.  Heterotrophic 
biomass  can  be  split  into  "nitrate  reducer"  and  "nitrite 
reducer" . The  nitrate  reducer  biomass  is  assumed  to  be 
capable  of  using  nitrate  oxygen  as  a terminal  electron 
acceptor.  This  division  of  heterotrophic  biomass  components 
need  not  be  interpreted  as  a suggestion  that  two  distinct 
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populations  of  denitrifiers  exist  (very  likely,  both  partial 
and  complete  denitrifiers  are  present) . Rather,  it  can  be 
interpreted  as  representing  the  respective  pools  of  the  enzyme 
complex,  i.e.,  nitrate  reductase  and  nitrite  reductase,  which 
are  distributed  among  the  heterotrophic  bacteria.  Autotrophic 
biomass  was  divided  into  ammonia  oxidizers  and  the  nitrite 
oxidizers.  The  ammonia  oxidizer  biomass  component  oxidizes 
ammonia  to  nitrite  and  nitrite  biomass  component  oxidizes 
nitrite  to  nitrate. 


Model  Matrix 

The  model  matrix  is  shown  in  Table  6-1.  The  components  in  the 
model  are  shown  at  the  top  of  the  matrix  and  are  described  on 
Table  6-2.  The  soluble  components  are  represented  by  the 
symbol  S,  whereas  particulate  components  are  represented  by 
symbol  X. 

Biological  processes  that  affect  the  components  in  the 
model  are  listed  in  the  far  left  column  of  the  table.  The 
corresponding  rate  expressions  are  given  in  the  far  right 
column.  Process  rates  are  denoted  by  the  symbol  p . An  index 
j is  assigned  to  each  of  the  processes  and  similarly  an  index 
i is  assigned  to  each  of  the  components. 

The  elements  in  the  matrix  are  the  stoichiometric 
coefficients,  v^.  They  signify  the  mass  relationships  between 
the  components  in  the  individual  processes. 


Table  6-1.  Process  kinetics  and  stoichiometry. 
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Table  6-2.  Components  of  the  model. 


i 

Symbol 

Description 

Units 

1 

xi 

Particulate  inert  substrate 

gCOD/m3 

2 

xs 

Particulate  substrate 

gCOD/m3 

3 

y 

"e, HN03 

Heterotrophic  biomass, 
nitrate  reducers* 

gCOD/m3 

4 

y 

d , HN02 

Heterotrophic  biomass, 
nitrite  reducers* 

gCOD/m3 

5 

y 

b, ANH3 

Autotrophic  biomass, 
ammonia  oxidizers* 

gCOD/m3 

6 

y 

B,AN02 

Autotrophic  biomass, 
nitrite  oxidizers* 

gCOD/m3 

7 

XND 

Particulate  biodegradable 
organic  nitrogen 

gN/m3 

8 

xp 

Particulate  products  from 
biomass  decay 

gCOD/m3 

9 

s> 

Soluble  inert  substrate 

gCOD/m3 

10 

Ss 

Soluble  substrate 

gCOD/m3 

11 

SN03 

Nitrate  nitrogen 

gN/m3 

12 

SN02 

Nitrite  nitrogen* 

gN/m3 

13 

SNH 

Ammonia  (ionized  + unionized) 
nitrogen 

gN/m3 

14 

SND 

Soluble  biodegradable  organic 

nitrogen  gN/m3 

15 

SALK 

Alkalinity  (HC03‘) 

moles/m3 

16 

S0 

Dissolved  oxygen 

gCOD/m3 

* Proposed  in  this  work. 
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The  sign  convention  used  in  the  matrix  is  positive  for 
production  and  negative  for  consumption.  The  concentrations 
of  all  organic  constituents,  whether  soluble  or  particulate, 
are  given  in  terms  of  Chemical  Oxygen  Demand  (COD) . Oxygen  is 
expressed  as  negative  COD.  The  concentration  of  biomass  and 
other  particulates  constituting  the  mixed  liquor  volatile 
suspended  solids  (MLVSS)  in  COD  units  should  be  multiplied  in 
dry  weight  units.  The  concentrations  of  nitrogen  forms  are 
given  in  terms  of  nitrogen. 

Some  of  the  stoichiometric  coefficients  in  the  proposed 
model  contain  conversion  factors  relating  COD  units  to 
nitrogen  units.  These  are  oxygen  utilization  by  autotrophs 
and  nitrate  utilization,  nitrite  production  and  utilization, 
and  alkalinity  production  by  heterotrophs  growing  under  anoxic 
conditions.  The  conversion  factors  can  be  derived  from  the 
half  reactions  for  nitrogen  conversion  and  the  half  reaction 
for  reduction  of  molecular  oxygen  to  water  (Grady  and  Lim, 
1980) : 


e‘ 

+ 

1/2 

no3‘ 

+ H+ 

1/2  N02'  + 

1/2 

h2o 

(6-1) 

e' 

+ 

1/3 

no2' 

+ 4/3  H+ 

= 1/6  N2  + 

2/3 

h2o 

(6-2) 

e* 

+ 

4/3 

H+ 

+ 1/6  N02' 

= 1/3  H20  + 

1/6 

nh4+ 

(6-3) 

e"  + 1/4  02  + H+  = 


1/2  H20 


(6-4) 
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The  derivation  of  each  conversion  factor  is  given  below: 

Oxygen  Utilization  by  Autotrophs  (coefficients  16,5  and  16.6) 
The  oxidation  of  ammonia  to  nitrate  occurs  in  two  steps 
in  the  proposed  model.  The  pertinent  reaction  for  the  first 
step  (oxidation  of  ammonia  to  nitrite)  is  obtained  by 
subtracting  equation  (6-4)  from  (6-3)  : 

1/6  NH4+  + 1/4  02  = 1/6  N02‘  + 1/6  H20  + 1/3H+  (6-5) 

In  this  reaction,  1 gram  NH4+-N  reacts  with  3.428  grams  02.  The 
pertinent  reaction  for  the  second  step  (oxidation  of  nitrite 
to  nitrate)  is  obtained  by  subtracting  equation  (6-4)  from  (6- 
1)  : 


1/2  N02'  + 1/4  02  = 1/2  N03'  (6-6) 

In  this  reaction,  1 gram  of  N02‘-N  reacts  with  1.143  grams  02. 
The  sum  of  these  two  coefficients  is  equal  to  the  single 
coefficient  of  4.57  that  is  used  in  the  IAWPRC  model  to  relate 
oxygen  consumption  to  the  growth  of  autotrophs. 

Nitrate  Utilization.  Nitrite  Production  and  Utilization,  and 
Alkalinity  Production  by  Heterotrophs  Growing  under  Anoxic 
Conditions  fcoef f icients  12.3;  11.3;  11.4;  15.4) 
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Nitrate  reduction  to  molecular  nitrogen  takes  place  in 
two  steps  in  the  proposed  model.  The  pertinent  reaction  for 
the  first  step  is  obtained  by  subtracting  equation  (6-4)  from 
(6-1) : 


1/2  NO3-  = 1/2  N02'  + 1/4  02  (6-7) 

This  reaction  shows  that  the  quantity  of  electrons  required  to 
reduce  1 gram  N03'-N  to  N02'-N  is  equivalent  to  1.143  grams  of 
oxygen.  The  pertinent  reaction  for  the  second  step  is 
obtained  by  subtracting  equation  (6-4)  from  (6-2) : 

1/3  N02‘  + 1/3  H+  = 1/6  N2  + 1/6  H20  + 1/4  02  (6-8) 

This  reaction  shows  that  the  quantity  of  electrons  required  to 
reduce  1 gram  N02'  - N to  N2  is  equivalent  to  1.714  grams  of 
oxygen.  The  sum  of  these  coefficients  is  equal  to  the 
coefficient  of  2.86  that  is  used  in  the  IAWPRC  model  to  relate 
nitrate  utilization  and  alkalinity  production  to  the  anoxic 
growth  of  heterotrophs . 

The  reaction  term  r;  for  a given  component  can  be  written 
by  adding  up  the  products  of  the  stoichiometric  coefficients 
v^.  and  process  rates  p.: 

ri  = E vij  Pj 

j= 0 


(6-9) 
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Model  Components 


The  model  given 

by  Henze 

et  al  (1986) 

contained 

13 

components , 

including 

two 

for 

biomass  (heterotrophic 

and 

autotrophic) 

and  two 

for 

inorganic  nitrogen 

(ammonia 

and 

nitrate) . The  model  proposed  herein  contains  16  components, 
including  4 for  biomass  and  3 for  inorganic  nitrogen  (Table  6- 
2) . The  biomass  and  inorganic  nitrogen  components  of  the 
proposed  model  are  described  next. 

Heterotrophic  biomass  is  divided  into  nitrate  reducing 
(XB  HN03)  and  nitrite  reducing  (XB  HN02)  components.  Both  grow  on 
readily  degradable  substrate,  using  oxygen  as  the  terminal 
electron  acceptor  under  aerobic  conditions  and  nitrate  or 
nitrite,  respectively,  as  the  terminal  electron  acceptors 
under  anoxic  conditions.  Growth  stops  under  anaerobic 
conditions,  in  which  neither  oxygen,  nitrate,  nor  nitrite  are 
present. 

Autotroph  bacteria  are  divided  into  ammonia  oxidizing 
(Xg  ANH4)  and  nitrite  oxidizing  (XB  AN02)  components.  Ammonia 
nitrogen  is  the  source  of  energy  for  the  growth  of  the  ammonia 
oxidizers,  as  well  as  source  of  nitrogen  for  biosynthesis  of 
both  heterotrophic  and  autotrophic  biomass.  Nitrite  nitrogen 
is  the  source  of  energy  for  the  growth  of  nitrite 
oxidizers.  Three  forms  of  inorganic  nitrogen  are 
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considered  in  the  proposed  model:  ammonia  (SNH4)  , Nitrite 

(sn02)  , and  nitrate  (SN03)  . 


Model  Processes 

The  proposed  model  considers  thirteen  processes,  in  comparison 
to  the  eight  described  by  Henze  et  al  (1990) . Those  differing 
from  the  original  model  are  covered  next. 

Aerobic  Growth  of  Nitrate  Reducing  Heterotrophs 

The  first  process  in  the  model  (row  1 in  Table  6-1)  is 
aerobic  growth  of  heterotrophic  biomass  capable  of  nitrate 
reduction.  Growth  occurs  at  the  expense  of  readily  degradable 
substrate  (Ss)  and  dissolved  oxygen  (S0)  , producing  X8HN03. 
Ammonia  nitrogen  (SNH4)  is  removed  from  solution  and 
incorporated  in  the  biomass.  Alkalinity  (SALK)  is  removed  in 
proportion  to  the  rate  of  ammonia  assimilation.  The  rate  of 
aerobic  growth  is  a function  of  the  concentrations  of  soluble 
substrate  and  dissolved  oxygen. 

Aerobic  Growth  of  Nitrite  Reducing  Heterotrophs 

The  aerobic  growth  of  the  nitrite  reducing  heterotrophs 
is  shown  in  row  2 of  Table  6-1.  This  process  produces  XB  HN02, 
but  is  other-wise  analogous  to  the  first  process. 
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Anoxic  Growth  of  Nitrate  Reducing  Heterotrophs 

The  anoxic  growth  of  heterotrophs  capable  of  nitrate 
reduction  is  considered  in  the  third  row  of  Table  6-1. 
Nitrate  nitrogen  (SN03)  acts  as  the  terminal  electron  acceptor 
in  this  environment.  Growth  occurs  at  the  expense  Ss  and  SN03, 
producing  XB  HM03  and  nitrite  (SNo2)  , SNH4  is  assimilated  for  cell 
biosynthesis  and  SALK  is  removed  in  proportion  to  the 
assimilation  of  SNH4.  A correction  factor,  r?GN03,  is  included 
in  the  process  rate  term  because  the  maximum  rate  of  substrate 
removal  is  found  to  be  less  under  anoxic  conditions  than  under 
aerobic  conditions.  Two  switching  functions  control  the 
relative  rate  of  growth  in  proportion  to  the  concentrations  of 
Ss  and  SN03.  A third  switching  function  turns  down  anoxic 
growth  in  proportion  to  the  concentration  of  SQ. 

Anoxic  Growth  of  Nitrite  Reducing  Heterotrophs 

The  anoxic  growth  of  heterotrophic  biomass  capable  of 
nitrite  reduction  is  considered  in  row  4.  This  process  is 
analogous  to  the  anoxic  growth  of  nitrate  reducing 
heterotrophs  except  that  SN02  acts  as  the  terminal  electron 
acceptor  in  the  anoxic  environment  and  XB  HN02  is  produced.  One 
additional  difference  is  that  SAL|C,  while  being  removed  in 
proportion  to  SNH4  assimilation,  is  also  produced  in  proportion 
to  the  growth  of  XB  HN02.  A correction  factor  (nGN02)  is  used  to 
account  for  the  slower  growth  using  nitrite  as  an  electron 
acceptor,  relative  to  growth  with  oxygen.  Switching  functions 
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for  Ss  and  SN02  turn  on  this  process  whereas  a switching 
function  for  So  turns  off  this  process. 


Growth  of  Ammonia  Oxidizers 

The  fifth  row  corresponds  to  the  aerobic  growth  of 
ammonia  oxidizers.  SNH4  is  oxidized  as  the  energy  source, 
consuming  SQ  and  producing  SN02  and  XB  ANH4.  Additional  SNH4  is 
assimilated  for  cell  biosynthesis.  SAL|C  is  consumed  for  cell 
biosynthesis  and  by  SNH4  oxidation.  Two  switching  functions 
control  the  relative  rate  of  growth  in  proportion  to  SNH4  and 

So‘ 

Growth  of  Nitrite  Oxidizers 

The  growth  of  nitrite  oxidizers  is  considered  in  row  6. 
SN02  is  oxidized  for  energy,  consuming  SQ  and  producing  SN03  and 
Xg  AN02.  SNH4  is  assimilated  and  SALK  is  consumed  for  cell 
biosynthesis.  Two  switching  functions  control  the  relative 
rate  of  growth  in  proportion  to  SN02  and  SQ. 

Decay  of  Nitrate  Reducing  Heterotrophs 

The  decay  of  nitrate  reducing  heterotrophic  bacteria  is 
considered  in  row  7.  The  rate  of  decay  is  dependent  on  the 
decay  coefficient,  bH  and  Xg  HNQ3.  Decay  results  in  the  loss 
of  XB  HN03  and  the  formation  of  particulate  products  (Xp)  , slowly 
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biodegradable  substrate  (Xs)  , and  particulate  degradable 
nitrogen  (XND)  . 

Decay  of  Heterotroohic  Nitrite  Denitrifiers 

The  decay  of  XB  HN02  is  considered  in  row  8.  This  process 
is  analogous  to  the  decay  of  nitrate  reducing  heterotrophs. 

Decay  of  Ammonia  Oxidizers 

The  decay  of  XB  ANH4  (row  9 of  the  matrix)  is  treated  in 
the  same  fashion  as  the  decay  of  heterotrophs.  The  decay 
coefficient,  bA,  is  different  from  that  of  the  heterotrophs. 

Decay  of  Nitrite  Oxidizers 

The  decay  of  the  nitrite  oxidizers  (row  10  of  matrix)  is 
analogous  to  the  decay  of  ammonia  oxidizing  bacteria 
(Hamilton,  1991)  . 

Ammonif ication  of  Soluble  Organic  Nitrogen 

Conversion  of  soluble,  degradable  organic  nitrogen  (SND) 
to  SNH4  is  considered  in  row  11.  This  process  is  first  order 
with  respect  to  SND  and  the  sum  of  XB  HN03  and  XB  HN02.  Protons 
are  consumed  on  an  equimolar  basis  with  SNH4  produced, 
increasing  the  alkalinity  of  the  system. 
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Hydrolysis  of  Entrapped  Organics 


Row 

12  of  the 

matrix  considers 

the 

hydrolysis 

of 

entrapped 

particulate 

degradable  (Xs) 

into 

Ss,  which 

is 

available  for  microbial  uptake.  The  rate  of  hydrolysis  is 
dependent  on  the  coefficient  KH  and  the  total  concentration  of 
heterotrophic  bacteria,  XB  HN02  and  XB  HN03.  A Contois-type 
switching  function  (Contois,  1959)  controls  the  relative  rate 
of  hydrolysis  in  proportion  to  two  ratios,  XS/XB  hno2  and 
XS/XB  HN03.  A second  function,  appearing  in  brackets  in  the  rate 
expression,  controls  the  relative  rate  of  hydrolysis  in 
proportion  to  electron  acceptor  concentrations.  A correction 
factor,  ?7h,  is  used  to  account  for  the  slower  rate  of 
hydrolysis  under  anoxic  conditions  as  compared  to  aerobic 
conditions . 

Hydrolysis  of  Entrapped  Organic  Nitrogen 

The  last  process  in  the  matrix  (row  13)  is  hydrolysis  of 
XND,  producing  SND.  The  rate  is  defined  by  the  hydrolysis 
reaction  for  Xs  and  the  ratio  of  XN0  to  Xs  in  the  reactor. 


Stoichiometric  Model  Parameters 
The  model  contains  seven  stoichiometric  parameters  (Table 
6-3).  Three  of  these  (fp,  iXB,  ixp)  are  unchanged  from  the 


original  IAWPRC  model . 
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Table  6-3.  Stoichiometric  parameters  for  model. 


Parameter 

Symbol 

Typical 

value 

Unit  note 

Heterotrophic  yield 
of  nitrate  reducers 

v 

iHN03 

0.67 

dimensionless 

(a) 

Heterotrophic  yield 
of  nitrite  reducers 

Y 

1 HN02 

0.67 

dimensionless 

(a) 

Autotrophic  yield 
of  ammonia  oxidizers 

Y 

ANH4 

0.212 

g COD/g  N 

(b) 

Autotrophic  yield 
of  nitrite  oxidizers 

Y 

iAN02 

0.029 

g COD/g  N 

(b) 

Fraction  of  biomass 
yielding  particulate 
products 

f 1 P 

0.08 

dimensionless 

(c) 

Mass  N/mass  COD 
in  biomass 

1XB 

0.086 

g N/g  COD 

(c) 

Mass  N/mass  COD  in 
products  from  biomass 
decay 

AX 

0.06 

g N/g  COD 

(c) 

(a)  Assumed  to  be  equal  to  the  value  for  overall 
heterotrophic  biomass  given  by  Henze  el  al.  (1986) 

(b)  Based  on  an  overall  yield  of  0.24  g cell  COD/g  NH4+-N 
oxidized  (Grady  and  Lim,  1980) 


(c)  Henze  et  al.  (1986) 
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The  autotrophic  yield,  YA , in  the  IAWPRC  model  was  the  yield 
for  the  overall  oxidation  of  ammonia  nitrogen  to  nitrate 
nitrogen.  The  autotrophic  yield  in  the  proposed  model  is 
divided  into  two  coefficients,  one  for  the  yield  of  ammonia 
oxidizers,  YANH4,  and  one  for  the  yield  of  nitrite  oxidizers, 
Y„uri, . The  values  for  these  coefficients  were  determined  from 
the  following  equations  (Grady  and  Lim,  1980) . Growth  of 
ammonia  oxidizers  (Nitrosomonas) : 

55NH4+  + 7602  + 109HC03‘  = 

C5H702N  + 57H20  + 104H2C03  + 54N02' 

(6-10) 


Growth  of  nitrite  oxidizers  (Nitrobacter) : 

400N02'+NH4++4H2C03+HC03'  + 19502  = C5H702N+3H20+4  00NC>3‘  (6-11) 


As  these  equations  show,  the  yield  of  Nitrosomonas  is 
0.212  g cell  COD  per  gram  NH4+-N  oxidized,  whereas  the  yield 
of  Nitrobacter  is  0.029  g cell  COD  per  gram  N02‘-N  oxidized. 
(Note:  1 gram  biomass  as  C5H702N  equals  1.42  grams  biomass  as 
COD).  The  sum  of  these  coefficients  is  0.24  grams  cell  COD 
per  gram  NH4+  - N,  as  used  in  the  IAWPRC  model.  The 
heterotrophic  yield  ( YH)  given  by  Henze  et  al,  (1986)  as 
typical  was  used  in  the  model  for  the  yield  of  nitrate 
reducing  heterotrophic  bacteria,  YHN03,  and  nitrite  reducing 
heterotrophic  bacteria,  YHN02. 
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Kinetic  Model  Parameters 

The  proposed  model  contains  22  kinetic  parameters,  as 
described  in  Table  6-4.  Default  values  for  these  parameters 
from  the  present  work  are  also  given  in  the  table.  The 
parameters  bA,  bH,  ka,  kh,  and  kx  are  unchanged  from  the 
IAWPRC  model.  The  other  parameters  are  discussed  next. 

Heterotrophic  bacteria  is  separated  into  two 
categories:  that  capable  of  nitrate  reduction  and  that 
capable  of  nitrite  reduction.  The  respective  maximum 
specific  growth  rates  are  symbolized  as  AHN03  and  AHN02*  T^e 
COD  half-saturation  coefficients  are  KSN03  for  nitrate- 
reducing  biomass  and  KSN02  for  nitrite-reducing  biomass. 

The  proposed  model  has  three  oxygen  switching  function 
parameters  that  control  heterotrophic  growth.  The 
coefficient  Kg  H is  the  oxygen  half-saturation  coefficient 
for  aerobic  growth  of  both  nitrate-reducers  and  nitrite- 
reducers.  The  coefficients  Kg  hno3  and  Kg  hno2  appear  in  the 
switching  functions  that  turn  down  the  anoxic  growth  of 
nitrate-reducers  and  nitrite-reducers,  respectively,  in 
proportion  to  SQ. 

There  are  two  nitrogen  half-saturation  coefficients  for 
the  anoxic  growth  of  heterotrophic  bacteria.  The 
coefficient  H appears  in  the  expression  for  anoxic  growth 
of  nitrate  reducers,  whereas  the  coefficient  H appears  in 
the  expression  for  anoxic  growth  of  nitrite  reducers. 
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Table  6-4.  Kinetic  parameters  in  the  proposed  model. 


Typical 


Parameter 

Symbol 

value 

Unit 

note 

Maximum  specific  growth 
rate  of  nitrate  reducers 

^HN03 

8.5 

d'1 

(a,d) 

Maximum  specific  growth 
rate  of  nitrite  reducers 

^HN02 

8.5 

d'1 

(a,d) 

COD  half  saturation 
coefficient  of 
nitrate  reducers 

^SNCtf 

20 

g COD/m3 

(b) 

COD  half  saturation 
coefficient  of 
nitrite  reducers 

KSN02 

5 

g COD/m3 

(c) 

Oxygen  switiching  function 
coefficient  for  anoxic  Kg  hno3 

growth  of  nitrate  reducers 

0.25 

-g  COD/m3 

(b) 

Oxygen  switiching  function 
coefficient  for  anoxic  Kg  hn02 

growth  of  nitrite  reducers 

0.20 

-g  COD/m3 

(c) 

Nitrate  half-saturation 
coefficient  for  nitrate 
reducers 

^03, H 

4.0 

g N03-N/m3 

(b) 

Nitrite  half-saturation 
coefficient  for  nitrite 
reducers 

^02,  H 

0.2 

g N02-N/m3 

(c) 

Decay  coefficient  for 
heterotrophic  biomass  bH 

(nitrate  & nitrite  reducers) 

1.09 

d'1 

(b,d) 
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Table  6-4.  Kinetic  parameters  in  the  proposed  model, 
(continued) 


Parameter 

Symbol 

Typical 

value  Unit 

Notes 

Correction  factor 
anoxic  growth  of  nitrate 
reducers 

^GN03 

0.8 

dimensionless 

(b) 

Correction  factor  for 
anoxic  growth  of  nitrite 
reducers 

^GW02 

0.8 

dimensionless 

(c) 

Correction  factor  for 
hydrolysis  under 
anoxic  conditions 

0.4 

dimensionless 

(b) 

Ammonif ication  Rate 

0.11 

m3/(g  celt  COO  * d) 

(b,  d) 

Maximum  specific 
hydrolysis  rate 

K 

5.2 

d'1 

(b,d) 

Half-saturation  Coefficient 
for  hydrolysis  of  slowly 
biodegradable  substrate 

0.05 

dimensionless 

(b,d) 

Maximum  specific  growth 
rate  of  ammonia  oxidizers 

^ANH4 

0.55 

d'1  (d,e,f) 

Maximum  specific  growth 
rate  of  nitrite  oxidizers 

^AN02 

0.72 

d'1  (d,  e,  f ) 

Ammonia  half -saturation 
coefficient  for  ammonia 
oxidizers 

0.58 

g NH4-N/m3 

(e,f) 

Nitrite-half  saturation 
coefficient  for  nitrite 
oxidizers 

^N02, A 

0.719 

g N02-N/m3 

(e,f) 

Oxygen  half-saturation 
coefficient  for  ammonia 

^0,NH4 

18.0 

-g  COD/m3 

(g) 

oxidizers 
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Table  6-4.  Kinetic  parameters  in  the  proposed  model, 
(continued) 


Parameter 

Typical 

Symbol  value 

Unit 

Notes 

Oxygen  half-saturation 
coefficient  for  nitrite 
oxidizers 

Kq.moz  96  -°  "<?  COD/m3 

(g) 

Decay  coefficient  for 
autotrophs  (ammonia  and 
nitrite  oxidizers) 

bA  0.1  d'1 

(b,d) 

(a)  Assumed  value, 

based  on  £H  from  Henze 

et  al . 

(1986) 

(b)  Henze  et  al.  (1986) 

(c)  Hamilton  (1991) 

(d)  Temp.  = 25  °C 

(e)  Experiments  of  this  research 

(f)  pH  = 7.0 

(g)  Assumed  value,  based  on  model  simulation 
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There  are  also  two  correction  factors  for  anoxic  growth  of 
heterotrophic  bacteria,  rjGN03  and  nGN02' 

Autotrophic  biomass  is  separated  into  ammonia  oxidizers 
and  nitrite  oxidizers,  each  with  its  own  maximum  specific 
growth  rate  (Aanh4  and  Aano2)  • There  are  two  nitrogen  half- 
saturation  coefficients  for  autotrophic  growth:  The 
coefficient  A appears  in  the  expression  for  growth  of 
ammonia  oxidizers,  whereas  the  coefficient  Kno2  A appears  in 
the  expression  for  growth  of  nitrite  oxidizers . There  are 
also  two  oxygen  half-saturation  coefficients  for  autotrophic 
growth,  K^,  and  K^. 

The  values  for  maximum  specific  growth  rates  and  half- 
saturation coefficients  for  ammonia  and  nitrite  oxidizers 
were  evaluated  in  chapter  3 and  chapter  4,  respectively. 

Two  oxygen  half-saturation  coefficients  for  autotrophic 
growth  were  determined  by  the  comparisons  between  simulation 
and  experiment. 


Stoichiometric  Coefficients 
Each  of  the  stoichiometric  coefficients  in  Table  6-1 
can  be  identified  by  its  position  in  the  matrix  as  mentioned 
earlier.  For  example,  the  coefficient,  1/YHN03,  has  the 
position  v101.  This  coefficient  gives  the  rate  of  change  in 
readily  degradable  substrate,  Ss,  when  multiplied  by  process 
rate  p.,.  The  origins  of  selected  stoichiometric 
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coefficients  (those  with  conversion  factors)  are  discussed 
next. 

Coefficient  v1f.  ; and  v16  6 

The  coefficients  v16  5 and  v16  6 give  the  effect  of 
autotrophic  growth  on  SQ.  The  value  3.428  in  the  first 
coefficient  represents  the  grains  02  (as  negative  COD) 
required  to  oxidize  1 gram  ammonia  nitrogen  to  nitrite 
nitrogen  in  the  absence  of  biosynthesis.  The  value  of  1.143 
in  the  second  coefficient  represents  the  quantity  of  02  (as 
negative  COD)  required  to  oxidize  1 gram  of  nitrite  nitrogen 
to  nitrate  nitrogen. 

Coefficient  v12  , and  v12  6 

The  coefficient  v12  3 and  v12  6 convert  the  anoxic  growth 
rate  of  nitrite  reducers  and  growth  rate  of  nitrite 
oxidizers  into  rates  of  nitrate  consumption  and  production, 
respectively.  Coefficient  v12  3 contains  the  value  1.143 
which  represents  the  quantity  of  electrons,  in  oxygen 
equivalents,  required  to  reduce  one  gram  of  nitrate  nitrogen 
to  nitrite  nitrogen. 

Coefficients  vi:  1 through  vi;  n 

The  coefficients  in  column  15  give  the  effect  of  growth 
and  ammonif ication  processes  on  alkalinity.  Each 
coefficient  contains  the  value  14  in  the  denominator,  which 
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converts  nitrogen  quantities  from  grams  to  moles.1  The 
value  1.714  in  coefficient  v15  4 represents  the  quantity  of 
electrons  (as  grams  02)  needed  to  reduce  1 gram  nitrite 
nitrogen  to  N2(g).  Finally,  coefficient  v15  5 contains  the 
value  7 in  the  denominator.  This  value  accounts  for  the  two 
moles  of  H+  produced  for  each  14  grams  ammonia  nitrogen 
oxidized. 

Coefficients  vn  3_and_y11  4 

Coefficients  vn  3 and  vn  4 convert  growth  rate  of  the 
nitrate  and  nitrite  reducers  to  rates  of  nitrite  production 
and  consumption,  respectively.  The  value  1.143  and  1.714 
were  discussed  previously. 

Computer  Simulation  Results 

Based  on  the  developed  process  model,  a computer 
program  was  written  using  QUICKBASIC.  Applying  kinetic  and 
stoichiometric  coefficients  attained  from  literature  reviews 
and  experiments  to  this  model,  computer  simulation  was 
carried  out  using  the  experiment  conditions  described 
earlier.  Since  the  influent  pump  was  operated  for  30 
minutes  each  hour,  a dynamic  model  using  the  fourth  order 
Runge-Kutta  method  was  developed.  In  these  simulations, 
IBM-compatible  computers  were  used.  Table  6-5  shows  steady- 


^he  second  term  in  v15  5 may  be  considered  to  have  the 
multiplier  2/14. 
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state  nitrogen  components  under  the  operating  conditions  of 
the  experiments  as  well  as  steady-state  experiment  data. 


Table  6-5.  The  experiment  and  simulation  results  under  same  operation  conditions  of 
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CHAPTER  7 

EVALUATION  OF  THE  PROPOSED  MODEL 
Introduction 

Two  different  simulations  were  carried  out  to  evaluate 
the  proposed  computer  model.  First,  a contact  stabilization 
system  was  simulated  using  the  effective  maximum  specific 
growth  rates  and  half-saturation  coefficients  for  Nitrosomonas 
and  Nitrobacter . estimated  from  fed-batch  experiments.  Two 
fed-batch  systems  were  simulated  to  verify  the  experiments. 

Statistical  Analysis  of  Experiment  and  Simulation  Data 
in  Contact  Stabilization  Process 

Alexander  et  al  (1980)  showed  a good  agreement  between 
experiment  and  simulation  data.  But  they  did  not  consider 
nitrite  nitrogen  as  a nitrogen  component  in  their  model. 

In  this  study,  percentages  of  oxidized  nitrogen  forms  in 
the  contact  tank  were  compared  to  display  the  agreement 
between  experiment  and  simulation  data.  The  percentage  of 
oxidized  nitrogen  forms  was  defined  as  the  ratio  of  the  sum  of 
oxidized  nitrogen  to  total  nitrogen  amount.  The  amount  of 
oxidized  nitrogen  is  represented  as  the  sum  of  nitrite  and 
nitrate  nitrogens,  whereas  the  sum  of  ammonia,  nitrite  and 
nitrate  nitrogens  is  described  as  the  total  nitrogen  amount. 
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Since  the  amount  of  organic  nitrogen  was  too  small  to  compare 
to  that  of  inorganic  nitrogen,  i.e.,  ammonia,  nitrite  and 
nitrate  nitrogen,  this  form  of  nitrogen  was  not  considered  in 
total  nitrogen  estimates.  The  percentage  of  nitrite  nitrogen 
was  expressed  as  the  ratio  of  the  amount  of  nitrite  nitrogen 
to  total  nitrogen  amount. 

These  percentages,  mentioned  previously,  were  not 
considered  in  reaeration  tank  since  only  nitrate  nitrogen  was 
in  reaeration  tank  due  to  a high  concentration  of  DO,  4 g/m3, 
and  longer  detention  time.  Though  a shorter  detention  time 
was  anticipated  at  higher  sludge  recycle  ratio,  the  detention 
time  for  the  reaeration  tank  was  almost  twice  as  much  as  that 
of  the  contact  tank. 

The  analysis  of  data  from  experiment  and  simulation  can 
be  carried  out  by  two  methods:  Find  a correlation  coefficient 
and  the  sum  of  the  difference  between  simulation  and 
experimental  data  under  same  operation  conditions.  The 
Execustat  (PWS-Kent  Publishing  Company,  1991)  was  used  to 
estimate  the  correlation  coefficients  among  percentages  of 
oxidized  nitrogen  and  nitrite  in  contact  tank. 

Figure  7-1  shows  oxidized  nitrogen  percentages  in  contact 
tank  at  different  sludge  retention  times.  The  percentage  of 
oxidized  nitrogen  increases  as  sludge  recycle  ratio  is  higher 
and  this  trend  is  demonstrated  regardless  of  SRTs. 
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Figure  7-1.  Percentages  of  oxidized  nitrogens  in  contact  tank 
according  to  SRTs  at  50%,  80%  and  110%  of  sludge 
recycle  ratios. 
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Figure  7-2.  Percentages  of  oxidized  nitrogens  in  contact  tank 
according  to  sludge  recycle  ratio  at  4 . 8 , 6.0  and 
8.0  days  of  SRTs. 
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The  correlation  coefficient,  r,  for  each  SRT  is  displayed  on 
the  corresponding  figure.  Although  r is  low,  0.7644,  at  8.0 
days  of  SRT,  experiment  and  simulation  data  show  a good 
correlation  between  recycle  ratio  and  percentage  of  oxidized 
nitrogens . 

Figure  7-2  shows  the  percentages  of  oxidized  nitrogen  at 
different  SRTs  and  presents  negative  correlation  coefficients 
between  experiment  and  simulation  data.  But  it  is  difficult 
to  conclude  that  simulation  and  experimental  data  have  a 
negative  relationship  since  the  maximum  difference  between  two 
values  is  less  than  15  percents. 

Figure  7-3  displays  the  relationships  between  nitrite 
nitrogen  percentages  in  contact  tank  and  sludge  recycle  ratios 
at  different  SRTs.  Experiment  and  simulation  data  show  the 
opposite  trend  at  4.8  days  of  SRT.  Nitrite  nitrogen 
percentages  in  contact  tank  are  generally  increased  with 
higher  sludge  recycle  ratio  at  SRTs  of  6.0  days  and  8.0  days 
though  they  have  poor  correlation  coefficients,  -0.1560  and 
0.2624,  respectively. 

The  relationships  between  SRT  and  nitrite  nitrogen 
percentage  in  contact  tank  were  determined  in  Figure  7-4.  As 
shown  on  these  figures,  the  correlation  coefficients,  r,  are 
all  negative  numbers.  The  maximum  absolute  difference  between 
two  data  is  13  percent  with  an  average  value  of  6.38%. 

Table  7-1  has  the  differences  between  simulation  and 
experiment  in  oxidized  nitrogen  percentage  under  same 
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Figure  7-3.  Percentages  of  nitrite  nitrogen  in  contact  tank 
according  to  SRTs  at  50%,  80%  and  110%  of  sludge 
recycle  ratios. 
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AT  50%  RECYCLE  RATIO 


r = -0.9888 


AT  80%  RECYCLE  RATIO 


r = -0.9840 


AT  1 10%  RECYCLE  RATIO 


r = -0.7896 
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4.5  5 5.5  6 6.5  7 7.5  8 8.5 

SRT,  DAY 

- EXPERIMENT  — SIMULATION 


Figure  7-4.  Percentages  of  nitrite  nitrogen  in  contact  tank 
according  to  sludge  recycle  ratio  at  4 . 8 , 6.0  and 
8.0  days  of  SRTs. 
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Table  7-1.  Percentage  differences  of  oxidized  nitrogens 

between  simulation  and  experiment  in  contact  tank. 


SRT 

20 

50 

80 

110 

140 

day 

% 

% 

% 

% 

% 

4.8 

11.42 

5.56 

3.07 

6 

1.33 

0.35 

1.23 

0.09 

1.33 

8 

0.30 

1.20 

0.12 

Average  = 2.36% 

Standard  Deviation  = 3.24% 


99 


operation  conditions.  The  difference  in  oxidized  nitrogen 
percentage  has  an  average  of  6.38%  with  a large  standard 
deviation,  5.65%.  This  large  standard  deviation  was  caused  by- 
two  operation  conditions  such  as  20%  and  50%  recycle  ratios  at 
6 days  of  SRT. 

An  average  in  difference  of  oxidized  nitrogens  was  2.36% 
as  shown  in  Table  7-2.  4.8  days  of  SRT  caused  a large 
standard  deviation,  3.24%. 

Simulation  of  Fed-Batch  Experiments 

Two  series  of  fed-batch  experiments  were  simulated  to 
evaluate  the  proposed  computer  model.  Exp.  #7/18A  and  B were 
simulated  as  shown  on  Figure  7-5.  Since  wasting  activated 
sludge  from  reaeration  tank  was  used  for  a source  of  biomass 
in  these  experiments,  the  initial  concentration  of 
Nitrosomonas  in  fed-batch  reactors  was  estimated  by  comparing 
the  MLSS  concentrations  between  the  reaeration  tank  and  the 
fed-batch  reactor.  On  the  date  of  the  experiment,  July 
18,1991,  the  reaeration  tank  had  1200  g/m3,  whereas  the  fed- 
batch  reactor  had  2 64  g/m3  of  MLSS  concentration.  The  initial 
Nitrosomonas  concentration  in  fed-batch  reactor  was  estimated 
to  be  78.4  g COD/m3  since  the  simulation  results  under  same 
operation  conditions,  i.e.,  8 days  of  SRT  with  110  % recycle 
ratio,  had  357.57  g COD/m3. 

Although  the  initial  Nitrosomonas  concentration  was 
estimated  as  78.4  g COD/m3,  the  simulation  in  fed-batch 
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Table  7-2.  Percentage  differences  of  nitrite  nitrogen 

between  simulation  and  experiment  in  contact  tank. 


SRT 

20 

50 

80 

110 

140 

day 

% 

% 

% 

% 

% 

4.8 

1.24 

1.72 

1.10 

6 

19.26 

11.87 

3.76 

7.94 

3.41 

8 

4.86 

2.19 

12.83 

Average  = 6.38% 

Standard  Deviation  = 5.65% 


101 


reactor  required  50.0  g COD/m3  as  its  initial  concentration  to 
simulate  experimental  data.  The  difference  of  two  initial 
concentrations  could  be  contributed  to  the  initial  conditions 
of  biomass  in  each  experiment.  Wasting  activated  sludge  from 
contact  stabilization  was  collected  for  the  period  of  one  day 
in  a bottle  in  refrigerator  and  this  refrigerated  sludge  was 
used  as  a source  of  biomass.  Refrigeration  may  have  decreased 
the  activity  of  microorganisms  and  therefore  caused  the 
difference  in  initial  biomass  concentration. 

In  the  case  of  Exp#  11/1A  and  B were  simulated  to 
evaluate  the  proposed  model.  On  November  1,  1991,  the  contact 
stabilization  system  was  operated  at  6 days  of  SRT  with  140% 
recycle  ratio.  Under  these  operation  conditions,  the  MLSS 
concentration  in  the  reaeration  tank  was  900  g/m3,  whereas  the 
fed-batch  reactor  had  125  g/m3.  Since  the  simulation  study 
under  the  same  operation  conditions  showed  291.6  g COD/m3  as 
the  concentration  of  Nitrobacter . its  coresponding 
concentration  in  the  fed-batch  reactor  was  estimated  as  40.5 
g COD/m3.  The  initial  concentration  of  Nitrobacter.  20.5  g 
COD/m3,  however,  showed  a good  correspondence  to  the 
experiment.  The  initial  acclimation  of  Nitrobacter  to  its  new 
environment  was  shown  in  the  experiment  and  simulation  as 
illustrated  on  Figure  7-6. 
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■ Experiment  — Simulation 


Figure  7-5.  experiment  and  simulation  data  in  study. 
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■ Experiment  — Simulation 


■ Experiment  — Simulation 


Figure  7-6.  Experiment  and  simulation  study  in  study. 


CHAPTER  8 
CONCLUSIONS 

The  following  conclusions  are  drawn  from  this  research: 

1.  Two  half-saturation  coefficients  for  ammonia  and 
nitrite  nitrogens  were  estimated  as  0.58  g NH4+-N/m3 
and  0.72  g N02'-N/m3,  respectively,  through  fed-batch 
experiments. 

2.  The  laboratory  fed-batch  experiments  gave  the 

effective  maximum  specific  growth  rates  for  ammonia 
and  nitrite  oxidizers  as  0.65  d'1  and  0.82  d'1, 

respectively . 

3.  Concentrations  of  near  30  g N02'-N/m3  in  reactor 

produced  substrate  inhibition  to  nitrite  oxidizer, 
Nitrobacter . in  fed-batch  experiments. 

4.  Sludge  Retention  Time  ( SRT)  did  not  show  much 

impact  on  nitrification  in  the  contact 

stabilization  process. 

5.  The  concentrations  of  nitrite  in  the  contact  tank 
showed  poor  correlations  according  to  either  SRT  or 
sludge  recycle  ratio. 

6.  Experiment  and  simulation  data  showed  a good 

correlations  for  the  percentage  of  oxidized 

nitrogens. 
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7.  There  was  poor  correlation  between  experimental  and 
simulation  data  for  the  percentage  of  nitrite 
nitrogen  in  the  contact  tank. 

8.  The  percentages  of  nitrite  nitrogen  in  contact  tank 
have  a poor  relationship  between  experiment  and 
simulation. 

9.  The  half-saturation  coefficients  study  showed  a 
good  correspondence  between  simulation  and 
experimental  data  despite  of  the  discrepancy  of  the 
initial  biomass  concentrations. 


APPENDIX  A 
EXPERIMENT  DATA 
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SRT  = 4.8  DAYS  WITH  110%  RECYCLE  RATIO 
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Rem  Program  for  Dynamic  Simulation  of  a WWTP  using  Runge-Kutta  method. 

Rem  Clarifier  is  assumed  to  be  perfect  separator 

Rem  Number  of  tanks  (in  series)  is  variable 

rem  Last  changes  : March  7,  1992 

rem  Define  Double  precision 

defdbl  a-h.o-z 

defint  i-n 

rem 

rem 

rem 

dim  dkl(4),o2con(4) 

dim  dkla(6),eqdo(6),phi(60,61),delta(60),codin(4) 
dim  codeff(4),codwaste(4),dsno(4) 

dim  XFeed(8,48),SFeed(8,48),x(8,4),s(8,4),dnue(16,13),rate(16,3),xx(50) 
dim  q(48),disso2(2) 

rem 

common  shared  yanh4.yano2,yhno3.yhno2,fp,dixb,dixp,_ 

dksno3,dksno2,dkohno3,dkohno2,dkno3h,dkno2h,dknh4a,dkno2a, 

dkonh4,dkono2,etah,etagno3,etagno2,dmuhno3,_ 

dmuhno2,bh,dka.dkh,dkx,ba,dmuanh4.dmuano2,dnue(),dkla(),_ 

eqdol,eqdo2,w,wl,w2,srt,x(),s(),  notanks,  r,_ 

disso2(),q(),eqdo(),vtotal,ihalf,XFeed(),SFeed(),deIt,_ 

q,ql,q2,q3,q4,qh,rate(),vl,v2 

rem 

rem  input.dat  contains  the  initial  values  at  time  t = 0 

rem  dyninp.dat  contains  time  varying  variables  ( value  at  0 and  48  th  hour 

rem  are  same) 

rem  output.dat  contains  daily  output 
rem  dynopt.dat  contains  dynamic  output 
rem 

open  "c:\qb87\model\sin_l.dat"  for  input  as  #1 
open  "c:\qb87\model\dIN48_l.dat"  for  input  as  #2 
open  "c:\qb87\model\outpt811.dat"  for  output  as  #5 
open  "c:\qb87\model\dynopt81.dat"  for  output  as  #6 
rem 

temp  = 25.0#  ’temp  is  temperature  in  degree  C. 

pH  = 7.0# 

rem 

rem  Stoichiometric  parameters  :- 
rem 

yanh4  = 0.212#  : yhno3  = 0.67#  : yhno2=0.67#  : fp  = 0.08#  : dixb  = 0.086# 

dixp  = 0.06#  : yano2  = 0.029# 

rem 

rem  Kinetic  parameters 

dksno3  = 20.0#  : dkohno3  = 0.25#  : dkno3h=4.0#  : dknh4a  = 0.58# 
dksno2=5.0#  : dkohno2  = 0.2#  : dkno2h  = 0.2#  : dkno2a  = 0.719# 
dkonh4  = 0.4# *45.0#  : dkono2  = 0.4#*240.0# 
etah  = 0.4#  : etagno3  = 0.7#  : etagno2  = 0.7# 
rem 

rem  Following  parameters  are  temperature  dependent 

dmuhno3  = 6.0#*exp(0.0693147#*(temp-20.0#)) 

rem  dmuhno3  = 5.2# 

dmuhno2  = 6.0#*exp(.069#  ’(temp- 15#)) 

rem  dmuhno2  = 5.2# 

bh  = 0.62#*exp(0. 1131402#  * (temp-20.0#)) 

rem  bh  = 0.5# 

dka  =0.08#*exp(0.0693147#*(temp-20.0#)) 
rem  dka  = 0.07 

dkh  =3.0#  *exp(0. 1098612#  * (temp-20.0#)) 
rem  dkh  = 2.4# 

dkx  =0.03#*exp(0.1098612#*(temp-20.0#)) 
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rera  dkx  = 0.024# 
rem  ba  =0.1# 

ba  = 0.09#*exp(0.1098612#*(temp-20.0#)) 

rem  dmuanh4  = 0.7#*exp(0.1283707#*(temp-20.0#)) 

dmuanh4  = 0 .551017# 

rem  dmuano2  = ((4.7*10~14)*exp((-9980)/(temp  + 273)))/(l  + (2.05*10~-9)/(10/'(-pH))  + (10~(-pH))/(1.66*10~-7)) 

dmuano2  = 0.720433# 

rem 

rem  X = > insoluble  constituents:- 

rem  1:  inert,  Xi  5:  autotrophic  biomass,  Xbanh4 

rem  2:  slowly  biodegradable,  Xs  (Ammonia-N  Oxidizing) 

rem  3:  heterotrophic  biomass,  Xbhno3  6:  autotrophic  biomass,  Xbano2 

rem  (Nitrate  Reducing)  (Nitrite  Dxidizing) 

rem  4:  heterotrophic  biomass,  Xbhno2  7:  Particulate  Product,  Xp 

rem  (Nitrite  Reducing)  8:  organic  nitrogen,  Xnd 

rem 

rem  S = > soluble  constituents:- 


rem 

1:  inert, 

Si 

5:  ammonia  nitrogen, 

Snh4 

rem 

2:  substrate, 

Ss 

6:  organic  nitrogen, 

Snd 

rem 

3:  nitrite-nitrogen, 

Sno2 

7:  alkalinity, 

Salk 

rem 

rem 

4:  nirtate-nitrogen, 

Sno3 

8:  dissolved  oxygen, 

So 

rem  i is  index  for  constituents 

rem  j is  index  for  process 

rem  k is  index  for  tanks  and  streams 

rem 

rem  Nues  are  stoichiometric  coeff.  in  matrix  in  Table  Nitrite  model 
for  i = 1 to  16 

for  j = 1 to  13 

dnue(i.j)  = 0.0# 

next 

next 

dnue(2,7)  = 1.0#-fp  : dnue(2,8)  = dnue(2,7)  : dnue(2,9)  = dnue(2,7) 
dnue(2,10)  = dnue(2,7)  : dnue(2,12)  = -1.0# 
dnue(3,l)  = 1.0#  : dnue(3,3)  = 1.0#  : dnue(3,7)  = -1.0# 
dnue(4,2)  = 1.0#  : dnue(4,4)  = 1.0#  : dnue(4,8)  = -1.0# 
dnue(5,5)  = 1.0#  : dnue(5,9)  = -1.0# 
dnue(6,6)  = 1.0#  : dnue(6,10)  = -1.0# 

dnue(7,7)  = fp  : dnue(7,8)  = fp  : dnue(7,9)  = fp  : dnue(7,10)  = fp 
dnue(8,7)  = dixb-fp*dixp  : dnue(8,8)  = dnue(8,7)  : dnue(8,9)  = dnue(8,7) 
dnue(8,10)=dnue(8,7)  : dnue(8,13)=-1.0# 

dnue(10,l)  = -1.0#/yhno3  : dnue(10,2)=-1.0#/yhno2  : dnue(10,3)  = -1.0#/yhno3 
dnue(10,4)  = -1.0#/yhno2  : dnue(10,12)  = 1.0# 
dnue(ll,3)  = (1.0#-yhno3)/(1.143#*yhno3)  : dnue(ll,5)  = 1.0#/yanh4 
dnue(ll,4)  = -(1.0#-yhno2)/(1.714#*yhno2)  : dnue(ll,6)  = -1.0#/yano2 
dnue(12,3)=-(1.0#-yhno3)/(1.143#*yhno3)  : dnue(12,6)  = 1.0#/yano2 
dnue(13,l)  = -dixb  : dnue(13,2)  = -dixb  : dnue(13,3)  = -dixb  : dnue(13,4)  = -dixb 
dnue(13,5)  = -dixb-1.0#/yanh4  : dnue(13,6)  = -dixb  : dnue(13,ll)  = 1.0# 
dnue(14,ll)  = -1.0#  : dnue(14,13)  = 1.0# 

dnue(15,l)  = -dixb/14.0#  : dnue(15,2)  = dnue(15,l)  : dnue(15,3)  = dnue(15,l) 
dnue(15,4)  = (1.0#-yhno2)/(14.0#*1.714#*yhno2)-dixb/14.0# 
dnue(15,5)  = -dixb/14.0#-1.0#/(7.0#*yanh4)  : dnue(15,6)  = -dixb/14.0# 
dnue(15,ll)  = 1.0#/14.0# 

dnue(16,l)  = -(1.0#-yhno3)/yhno3  : dnue(16,2)  = -(1.0#-yhno2)/yhno2 
dnue(16,5)  = -(3.428#-yanh4)/yanh4  : dnue(16,6)  = -(1.143#-yano2)/yano2 


REM 

Rem  Define  number  of  tanks  in  system  (including  anoxic  basin) 
NoTanks  = 2 

rem  define  klal,kla2  (mass  transfer  coeff.  for  02  ) here 
dklal  = 0.0#  : dkla2  = 102.0#  ’per  day 
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rem 

rem  define  saturation  DO  concentration 
eqdol  = 7.56#  : eqdo2  = 8.637#  ’gm  02/m3 

rem 

disso2(l)  = 4.0#  : disso2(2)=4.0#  'Dissolved  02  cone,  in  aeration  basin 
rem 

rem  define  tank  volumes  here 

vl  =5.0#  : v2  = 2-5# 

rem 

w=0.6# 
r=l.l# 
q = 12.5# 

srt  = 8.0#  ’ Solid  Retention  Time  in  days 

rem 

rem  read  initial  values  at  t = 0 
rem 

for  ii  = 1 to  8 

input#l,XFeed(ii,0),x(ii,l)>x(ii,2) 

next  ii 

for  ii  = 1 to  8 

input#l,SFeed(ii,0),s(ii,l),s(ii,2) 

next  ii 
rem 

rem  read  dynamically  varying  variables 
rem 

for  j = 0 to  24 

input#2,ihalf,q(j),XFeed(l,j),XFeed(2,j),XFeed(3,j),XFeed(4,j),XFeed(5  j),XFeed(6,j),XFeed(7,j),XFeed(8,j), 
SFeed(lj),SFeed(2j),SFeed(3j),SFeed(4,j),SFeed(5j),SFeed(6,j),SFeed(7,j),SFeed(8j) 

next  j 
rem 

print#5,  " SRT  = ";srt;"  days" 

print#5,  " Initial  values 

ihalf =0 

imin  = 0 

call  Printlnput 

rem  tank  1&2  constituent  concentration  are  mapped  into  array  XX 
rem  x(k,l)  & x(k,2),k  = 1..6  go  in  xx(l)..xx(6)  & xx(14)..xx(19)  respectively, 
for  i = 1 to  8 

xx(i)  = x(i,l) 
xx(i  + 16)  = x(i,2) 

next  i 

rem  s(k,l)  & s(k,2),k=1..7  go  in  xx(7)..xx(13)  & xx(20)..xx(26)  respectively, 
for  i = 1 to  8 

xx(i  + 8)  = s(i,l) 
xx(i  + 24)  = s(i,2) 

next  i 
ncrvar=32 
t = 0.0# 

delt  = 1.0#/24.0#/60.0#  'days 
tmax  = 20.01#  'days 
icode  = 1 

call  PrintParameters(icode) 

call  RungeKutta(novar,t,delMn>ax,xx()) 

close#l 

close#2 

close#5 

close#6 

end 


sub  Printlnput  static 


122 


print#5,"  " 

print#5, 'COMPONENT, "INPUT","  TANK  1","  TANK  2" 

PRINT#5, ' 

PRINT#5, " Xi  ”,XFeed(l,ihalf),x(l,l),x(l,2) 

PRINT#5,  * Xs  ",XFeed(2,ihalf),x(2,l),x(2,2) 
print#5, " Xbhno3  ",XFeed(3,ihalf),x(3,l),x(3,2) 
print#5,  " Xbhno2”,XFeed(4,ihalf),x(4,l),x(4,2) 
print#5,  " Xbanh4  ",XFeed(5,ihalf),x(5,l),x(5,2) 
print#5,  ' Xbano2",XFeed(6,ihalf),x(6,l),x(6,2) 
print#5,  ’ Xp  ",XFeed(7,ihalf),x(7,l),x(7,2) 
print#5, ' Xnd  ",XFeed(8,ihalf),x(8,l),x(8,2) 
print#5, " Si  ",SFeed(l,ihalf),s(l,l),s(l,2) 
print#5,  ’ Ss  ",SFeed(2,ihalf),s(2,l),s(2,2) 
print#5,  ’ Sno2  ",SFeed(3,ihalf),s(3,l),s(3,2) 
print#5, " Sno3  ",SFccd(4,ihalf),s(4,l),s(4,2) 
print#5,  ' Snh4  ",SFeed(5,ihalf),s(5,l),s(5,2) 
print#5,  ' Snd  ",SFeed(6,ihalf),s(6,l),s(6,2) 
print#5, ' Salk",SFeed(7,ihalf),s(7,l),s(7,2) 
print#5, " So  ”,SFeed(8,ihalf),s(8,l),s(8,2) 
print#5, ' " 

print#5, ' " 

end  sub 


sub  EvaluateFunctions(xx(l),f(l))  static 
rem  dim  rate(16,3) 

rem  map  composition  array  XX  back  into  soluble  & particulate  arrays, 
ql  = r*q-w 
q2  = qh  + r*q-w 
q4  = r*q 


for  k = 1 to  8 

x(k,l)  = xx(k) 
x(k2)  = xx(k+ 16) 

next  k 

for  k = 1 to  8 

s(k,l)  = xx(k  + 8) 
s(k^)  = xx(k  + 24) 

next  k 

for  i = 1 to  8 

x(i,4)  = q2/q4*x(i,2) 

next  i 

for  i = 1 to  8 

s(i,4)  = s(i,2) 

next  i 

call  Determine  Rates(notanks,x(),s(),rate()) 
for  i = 1 to  8 

f(i)  = ((q4*x(i,4)-ql*x(i,l)-w*x(i,l))  + rate(i,l)*vl)/250.0# 
f(i  + 8)  = ((q4*s(i,4)-ql*s(i,l)-w*s(i,l))  + rate(i  + 8,l)*vl)/250.0# 
f(i  + 16)  = ((qh  * XFeed(i,ihalf)  + ql*x(i,l)-q2*x(i,2))  + rate(i,2)*v2)/250.0# 
f(i  + 24)  = ((qh*SFeed(i,ihalf)  + ql*s(i,l)-q2's(i,2))  + rate(i  + 8,2)*v2)/250.0# 
rem  print, t,ql,q2,q4 
next  i 

f(16)  = s(8,l)-disso2(l) 

rem  f(15)  = (q(ihalf)*SFeed(7,ihalf)-ql*s(7,l)  + rl*s(7,2)  + r2*s(7,4))/vl  + dklal*(eqdol-s(7,l)) 
rem  +rate(15,l) 

REM  f(31)  = (ql*s(7,l)-q2*s(7,2)-rl*s(7,2))/v2  + dkla2*(eqdo2-s(7,2))_ 

REM  + rate(15,2) 

f(32)  = s(8,2)-disso2(2)  ’ set  DO  in  aeration  basin  at  disso2 


end  sub 
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sub  DetermineRates(notanks.x(2),s(2),rate(2))  static 
dim  switchl(3),switch2(3),switch3(3),switch4(3),switch5(3),switch6(3),ro(13,3) 
dim  switch7(3),switch8(3),switch9(3),switchl0(3),switchll(3),switchl2(3) 
for  k = 1 to  notanks 
switchl(k)  = s(2,k)/(dksno3  + s(2,k)) 
switch2(k)  = s(8,k)/(dkohno3  + s(8,k)) 
switch3(k)  = s(2,k) /(dksno2  + s(2,k)) 
switch4(k)  = s(8,k)/(dkohno2  + s(2,k)) 
switch5(k)  = dkohno3/(dkohno3  + s(8,k)) 
switch6(k)  = s(4,k)/(dkno3h  + s(4,k)) 
switch7(k)  = dkohno2/(dkohno2  + s(8,k)) 
switch8(k)  = s(3,k)/(dkno2h  + s(3,k)) 
switch9(k)  = s(5,k)/(dknh4a  + s(5,k)) 
switchlO(k)  = s(8,k)/(dkonh4  + s(8,k)) 
switchll(k)  = s(3,k) /(dkno2a  + s(3,k)) 
switchl2(k)  = s(8,k) /(dkono2  + s(8,k)) 


ro(l,k)  = dmuhno3*switchl(k)*switch2(k)*x(3,k) 

ro(2,k)  = dmuhno2*switch3(k)*switch4(k)*x(4,k) 

ro(3,k)  = dmhno3*switchl(k)*switch5(k)*switch6(k)*etagno3*x(3,k) 

ro(4,k)  = dmuhno2*switch3(k)*switch7(k)*switch8(k)*etagno2*x(4,k) 

ro(5,k)  = dmuanh4*switch9(k)*switchl0(k)*x(5,k) 

ro(6,k)  = dmuano2*switchll(k)*switchl2(k)*x(6,k) 

ro(7,k)  = bh*x(3,k) 

ro(8,k)  = bh*x(4,k) 

ro(9,k)  = ba*x(5,k) 

ro(10,k)  = ba*x(6,k) 

ro(ll,k)  = dka*s(6,k)*x(3,k)  + dka*s(6,k)*x(4,k) 

ro(12,k)  = dkh*x(2,k)/((x(3,k)  + x(4,k))*dkx  + x(2,k))*(switch2(k)  + etah*switch5(k) 
*switch6(k)  + etah*switch7(k)*switch8(k))*(x(3,k)  + x(4,k)) 
ro(13,k)  = dkh*x(8,k)/((x(3,k)  + x(4,k))*dkx  + x(2,k))*(switch2(k)  + etah*switch5(k) 
*switch6(k)  + etah*switch7(k)*switch8(k))*(x(3,k)  + x(4,k)) 


for  i = 1 to  16 

rate(i,k)  = 0.0# 
for  j = 1 to  13 

rate(i,k)  = rate(i,k)  + dnue(i,j)*ro(j,k) 

next  j 

next  i 
next  k 
end  sub 

sub  RungeKutta(n,t,delt,tmax,xx(l))  static 

dim  delx(60),xxl(60),xx2(60),xx3(60),dkl(60),dk2(60),dk3(60),dk4(60) 
a = 1.0#:b  = 2.0#:c  = 2.0#:d  = 1.0# 
index  = 0 

jindex  = 0 ’jindex  accounts  day 
imin  = 0 

while  (t  < = tmax) 

rem  ihalf =int((t-int(t))*48.0#  + 0.001#) 

imin  = int((t-int(t))  *24.0#  *60.0#  + 0.001#) 

Iflow= int((imin  + 30)/30) 
flowstat  = Iflow  mod  2% 
if  flowstat  = 0#  then 
qh  = 0.0# 
else 

qh  = 25.0# 
end  if 

call  EvaluateFunctions(xx(),dkl()) 
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for  j = 1 to  n 

xxl(j)  = xx(j)  + delt*dkl(j) 

next  j 

call  EvaluateFunctions(xxl(),dk2()) 
for  j = 1 to  n 

xx2(j)  = xx(j)  + delt*dk2(j) 

next  j 

call  EvaluateFunctions(xx2(),dk3()) 
for  j = 1 to  n 

xx3(j)  = xx(j)  + delt  *dk3(j) 

next  j 

call  EvaluateFunctions(xx3(),dk4()) 
for  j = 1 to  n 

delx(j)  = delt*(a*dkl(j)  + b*dk2(j)  + c*dk3Q)  + d*dk4(j))/6.0# 
xx(j)  = xx(j)  + delx(j) 

next  j 

rem  map  composition  array  XX  back  into  soluble  & particulate  arrays, 
for  k = 1 to  8 

x(k,l)  = xx(k) 
x(k,2)  = xx(k  + 16) 

next  k 

for  k = 1 to  8 

s(k,l)  = xx(k  + 8) 
s(k,2)  = xx(k  + 24) 

next  k 

for  i=  1 to  8 

x(i,4)  = q2/q4*x(i,2) 
s(i,4)  = s(i,2) 

next  i 

rem  print,t,s(3,2),s(4,2),s(5,2) 
print,  jindex 

print  using  "####.##";  t, 
print  using  "####.##";  qh, 
print  using  "#####.#";  s(3,l), 
print  using  "#####.#";  s(3,2), 

print  using  "#####.#";  s(4,l), 
print  using  "#####.#";  s(4,2), 

print  using  "#####.#";  s(5,l), 
print  using  "#####.#";  s(5,2) 

rem  following  code  provides  dynamic  output  4 times  a day  (with  ifreql  = 600) 
rem  for  more  frequent  output  decrease  ifreql  and  vice  versa 
rem  if  (index  > 725)  and  (index  < 735)  then 
if  (index  = 730)  then 

rem  if  t- int(t  + . 0001#) <.0003#  then 

rem  if  (int(t*24.0#*60.0#  + 0.0001#)  mod  ifreql)  = 0 then 
rem  map  composition  array  XX  back  into  soluble  & particulate  arrays, 
for  k = 1 to  8 

x(k,l)  = xx(k) 
x(k,2)  = xx(k  + 16) 

next  k 

for  k = 1 to  8 

s(k,l)  = xx(k  + 8) 
s(k,2)  = xx(k  + 24) 

next  k 

for  i = 1 to  8 

x(i,4)  = q2/q4*x(i,2) 
s(i,4)  = s(i,2) 


next  i 
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rem  pnnt,t,s(3,2),s(4,2),s(5,2) 

print  using  '####.##";  t, 
print  using  "####.##";  qh, 
print  using  "#####.#”;  s(3,l), 
print  using  "#####.#";  s(3,2), 

print  using  "#####.#";  s(4,l), 
print  using  "#####.#";  s(4,2), 

print  using  "#####.#";  s(5,l), 
print  using  "#####.#";  s(5,2) 

print#6,t,s(3,l),s(4,l),s(5,l),s(3,2),s(4,2),s(5,2) 

rem  call  PrintOutput(t,s(),x()) 
rem  end  if 
end  if 

rem  Following  code  prints  daily  output  (with  ifreq2  = l) 
rem  for  less  frequent  output, say  every  3 days,  make  ifreq2  = 3 
ifreq2  = 19 

rem  ifq2  = 1400*ifreq2 

if  (jindex  > 0 ) then 
if  (jindex  mod  ifreq2)  = 0 then 
if  (index  = 1439)  then 
rem  ifreq2  = 2400.0# 

rem  if  t- int(t  + .0001#)  < .0003#  then 
rem  if  (int(t*2400  + 0.00001#)  mod  ifreq2)  = 0 then 
icode  = 2 

call  PrintParameters(icode) 
print#5,"  " 
print#5,"  " 

rem  print#5,"  Day  ";int(t  + 0.0001#) 

print#5,"  Day  ";jindex,  "SRT  = ";SRT,  "r  =";r 
call  Printlnput 
end  if 
end  if 
end  if 

t = t + delt 
index  = index  + 1 
if  (index  =1440)  then 
index  = 0 

jindex = jindex + 1 


rem 

if  (jindex  = 21)  then 

rem 

jindex  = 0 

rem 

end  if 

end  if 

rem 

if  (imin  = 1440)  then 

rem 

imin  = 0 

rem 

end  if 

wend 

end  sub 

rem  sub  PrintOutput(t,s(2),x(2))  static 

rem  print#6,t;s(3,l);s(4,l);s(5,l);s(3,2);s(4,2);s(5,2) 

rem  end  sub 


sub  PrintParameters(icode)  static 
if  (icode  = 1)  then 

print#5,  " Kinetic  and  Stoichiometric  Parameters  for  Heterophic  Biomass'' 
print#5,  " 

print#5,  " Parameter  Nitrite  Nitrate  Units" 

print#5,  " reducer  reducer  " 
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print#5,  " 
print#5,  " 
print#5,  " 
print#5,  " 
print#5,  " 
print#5,  " 
print#5,  " 
print#5,  " 
print#5,  " 
print#5,  " 
print#5,  " 
print#5,  " 
print#5,  ” 
print#5,  " 
print#5,  " 
print#5,  " 
print#5,  " 
print#5,  " 
print#5,  " 
print#5,  " 
print#5,  " 
print#5,  " 
print#5,  " 
print#5,  " 
print#5,  " 
print#5,  " 
print#5,  " 
print#5,  " 
print#5," 
print#5, " 
print#5,  * 
print#5,  " 
print#5,  " 
print#5,  " 
print#5,  " 
print#5,  " 
print#5,  " 
print#5,  " 
print#5,  " 
print#5,  " 
print#5,  " 
print#5,  " 
print#5,  " 


",dmuhno2," 


",dksno2, 


",dkohno2," 


",yhno2," 


",bh," 


1/day" 

",dmuhno3,"l/day" 

gm  COD/m3" 
",dksno3,"  gm  COD/m3 
gm  02/m3" 
,dkohno3,"gm  02/m3 
gm/gm" 
gm/gm" 

1/day" 

1/day" 


",yhno3," 


",bh, 


",etagno3 


',dkno3h,’ 


gm  N/m3  ' 
gm  N/m3  ' 


Mu  max 
Mu  max 
Ks  COD 
Ks  COD 
Ko  02 
Ko  02 
Yield 
Yield 
b decay 
b decay 

anoxic  growth  ",etagno2 
factor 

Ks  N02  ",dkno2h," 

Ks  N03 
Hydrolysis  rate", dkh,’ "1/day  " 

Hydrolys.  satur.",dkx,"gm  COD/gm  COD" 
ratio 

Anoxic  hydrolys.", etah 
factor  " 

Ammonific.  rate.",dka,"  m3/gmCOD/day" 
fraction  panic.", fp,"  gm  COD/gm  COD" 
produced  " 

N in  biomass  “,dixb,"  gm  N/gm  COD" 

N in  particulate", dixp,"  gm  N/gm  COD 

product  " 


Kinetic  and  Stoichiometric  Parameters  for  Autotrophic  Biomass" 

Parameter  Ammonia  Nitrite  Units" 

Oxidizer  Oxidizer  " 

N 

",dmuanh4,dmuano2, "1/day" 

",dknh4a,"  gm  N/m3  " 

",dkno2a,"gm  N/m3  " 

",dkonh4,"  gm  N/m3  " 

",dkono2,"gm  N/m3  " 

",yanh4,"  gm/gm" 

",yano2,"  gm/gm" 

",ba,"  1/day" 

",ba,"  1/day" 


Mu  max 
Ks  NH4 
Ks  N02 
Ko  NH4 
Ko  N02 
Yield 
Yield 
b decay 
b decay 


else 

print#5,  " ’ 

print#5,  " " 

print#5,  " Process  Description  " 

print#5,  " " 

print#5,  " Number  of  Reactors  = ",2 
print#5,  " Solid  Retention  Time  = ",srt,"  days" 
print#5,  " Diss.  02  in  Reaeation  Tank=  ",disso2(l),"g/m3" 
print#5,  " Diss.  02  in  contact  tank  = ",disso2(2),"g/m3" 
print#5,  " Average  Flow  Rate  = ",q,"m3/day" 
print#5,  " Activ.  sludge  recyc.  ratio  = ",r 
print#5,  " Volume  of  reactor  1 = ",vl,"m3" 

print#5,  " Volume  of  reactor  2 = ",v2,"m3" 

print#5, " " 
end  if 
end  sub 


APPENDIX  C 

EXPERIMENT  AND  SIMULATION  GRAPHS 
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AT  SRT  = 4.8  days 
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AT  50%  RECYCLE  RATIO 


AT  80%  RECYCLE  RATIO 


i 1 1 1 r 


AT  110%  RECYCLE  RATIO 
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